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ISSUE:  Corps  dredging  operation  elements 
must  often  assess  the  retention  of  material  at 
dredged  material  disposal  sites  to  determine 
site  capacity  for  future  dredging  work.  This  is 
routinely  done  by  bathymetric  surveys  of  the 
site.  In  some  cases,  predictions  must  be  made 
of  sediment  erosion  processes  that  might 
occur  at  the  site  under  certain  assumed  condi¬ 
tions,  such  as  for  a  different  sediment  type. 
Erosion  depends  strongly  on  bed  shear  stress. 
A  variety  of  approaches  are  available  to  calcu¬ 
late  boundary  layer  properties  such  as  shear 
stress.  The  issue  is  to  determine  which 
method  is  most  appropriate  when  conditions 
are  complicated  by  tidal  influences. 

RESEARCH:  The  Dredging  Research  Pro¬ 
gram  (DRP)  includes  guidance  and  proce¬ 
dures  for  predicting  movements  of  dredged 
material  placed  in  open  water.  An  important 
part  of  DRP  research  is  to  check  predictive 
techniques  against  data,  as  was  done  here. 


SUMMARY :  Conditions  that  can  be  en¬ 
countered  in  estuarine  areas  include  tidally  re¬ 
versing  flows,  density  stratification,  sus¬ 
pended  sediment  contributions,  and  combined 
current  and  wind-wave  effects.  Possible  ap¬ 
proaches  to  making  boundary  layer  calcula¬ 
tions  include  modified  steady-flow  equations, 
variable  eddy  viscosity  and  mixing  length 
equations,  and  turbulence  transport  equations. 
Interim  recommendations  for  calculation 
methods  are  based  on  the  situation  and  the 
computational  power  available. 

AVAILABILITY  OF  REPORT:  The  report 
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FOREWORD 


In  1932  Sir  Horace  Lamb  said,  "I  am  an  old  man  now,  and  when  I  die  and 
go  to  heaven  there  are  two  matters  on  which  I  hope  for  enlightenment.  One  is 
quantum  electrodynamics,  and  the  other  is  turbulent  motion  of  fluids.  About 
the  former  I  am  really  rather  optimistic."*  Fifty-eight  years  later,  we  have 
little  more  reason  for  optimism  than  did  Sir  Horace.  We  have  made  progress  in 
the  sense  that  there  are  several  more  options  for  mathematically  describing 
fluid  turbulence  and  its  offspring,  fluid  stresses;  however,  these  options  are 
less  the  result  of  improved  physical  understanding  than  of  hugely  increased 
computational  capabilities. 

We  share  Sir  Horace's  hope  for  enlightenment,  but  confess  that  the  tech¬ 
niques  recommended  here  do  not  provide  that  enlightenment.  They  simply  pro¬ 
vide  means  of  computing  stresses  that  are  the  best  currently  available 
engineering  approximations. 


*  R.  J.  Donnelly.  1988  (Nov).  "Superfluid  Turbulence,"  Scientific  American 

p  100. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply 
cubic  feet 
degrees  (angle) 
feet 


- E* _ 

0.02831685 

0.01745329 

0.3048 


To  Obtain 
cubic  metres 
radians 
metres 
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SUMMARY 


Calculation  of  sediment  entrainment  and  transport  rates  at  open  water 
dredged  material  placement  sites  in  tidal  waters  requires  accurate  descrip¬ 
tions  of  velocity  profiles  and  stresses  on  the  sediment  bed.  Standard  equa¬ 
tions  for  calculations  in  steady,  uniform  flows  can  give  misleading  results 
when  applied  in  tidal  waters.  For  example,  standard  log  profiles  adapted  from 
steady- flow  equations  can  generate  bed  shear  stress  values  that  are  wrong  by  a 
factor  of  two  or  more  and  are  in  the  wrong  direction  as  well. 

Complicating  factors  in  tidal  water  environments  include  tidally  revers¬ 
ing  flows  in  which  accelerations  and  decelerations  strongly  affect  the  veloc¬ 
ity  profile  and  bed  shear.  Density  stratification,  common  in  estuaries,  can 
likewise  generate  profiles  that  drastically  differ  from  those  seen  in  homo¬ 
geneous  flows.  Suspended  sediment  contributions  and  deformable  mud  beds  also 
interact  with  the  flow,  altering  the  stresses  experienced  by  sediment  beds. 
Combined  current  and  wind  wave  effects  lead  to  stresses  on  the  bed  that  are 
not  simply  linear  combinations  of  their  individual  effects. 

Possible  approaches  to  making  calculations  in  this  environment  include 
(a)  modified  applications  of  the  steady-flow  equations,  (b)  time-  and  space- 
varying  eddy  viscosity  ana  mixing  length  equations,  and  (c)  turbulence 
transport  equations . 

Interim  recommendations  for  calculation  methods  are  based  on  the  situa¬ 
tion  and  the  computational  power  available: 

a.  For  simple  calculations,  equations  using  a  spatially  variable  eddy 
viscosity  are  recommended.  Where  density  stratification  and/or 
sediment  suspension  conditions  are  significant,  modifications  to 
the  eddy  viscosity  equations  are  justified.  Different  sets  of 
equations  are  recommended  for  tidal  currents  alone  and  for  tidal 
currents  with  short-period  waves  present. 

b.  For  three-dimensional  numerical  modeling,  a  turbulence  transport 
model  using  a  two-equation  closure  is  recommended. 

Laboratory  and  field  tests  are  needed  to  confirm  these  recommendations 
and  to  extend  them  to  more  realistic  environmental  conditions. 
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SOmDARY-STRESg£S,.A^£.  VELOCITY  PROFILES  IN  ESTUARINE  FLOWS 
INTERIM  CALCULATION  METHODS 


PART  I :  INTRODUCTION 

Objectives 

1.  The  primary  objective  of  this  work  is  to  provide  improved  descrip¬ 
tions  of  boundary  layer  velocity  profiles  and  stresses  in  estuarine  flows  such 
that  more  accurate  and  reliable  computations  can  be  made  for  dredged  material 
erosion  and  transport.  The  purpose  of  this  report  is  to  establish  the  best 
methods  currently  available  to  calculate  estuarine  velocity  profiles  and 
stresses.  A  subsequent  report  will  provide  updated  recommendations. 

Background 

2.  Accurate  computation  of  sediment  deposition,  erosion,  entrainment, 
and  transport  requires  realistic  calculation  of  near-bed  flows  and  effective 
stresses.  Deposition  and  erosion  rate  expressions  for  cohesive  sediments 
typical  of  estuaries  include  the  effective  shear  stress  exerted  on  the  bed  as 
a  basic  variable.  The  near-bed  velocity  profile  and  turbulent  exchanges  con¬ 
trol  the  rate  of  transport  of  sediments  traveling  close  to  the  bed. 

3.  While  this  important  subject  has  been  thoroughly  explored  for 
steady,  constant -density  flows,  it  has  been  largely  neglected  for  estuarine 
flows,  which  are  typically  unsteady  (due  to  astronomical  and  wind  tides, 
surges,  and  short-period  waves),  nonuniform,  and  variable  density  (due  to 
freshwater -saltwater  interaction  plus  temperature  and  suspended  sediment). 
Near-bed  processes  for  short-period  waves  have  been  examined  somewhat  more 
fully  because  of  the  obvious  inadequacies  of  traditional  approaches,  but  the 
mathematical  descriptions  of  those  are  still  Inadequate  to  use  in  computations 
involving  estuaries.  Estuarine  flows  have  usually  been  described  using  the 
traditional  steady-flow  methods  without  modification,  despite  the  fact  that 
they  can  give  erroneous  or  even  nonsensical  results. 

4.  As  part  of  an  earlier  attempt  to  develop  computational  procedures 
for  a  two-dimensional  numerical  model  of  sediment  transport,  some  ad  hoc 
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solutions  to  this  problem  were  developed  (Thomas  and  McAnaliy  1985) .  The  work 
described  here  is  an  attempt  to  improve  upon  that  approach  and  sec  a  framework 
for  future  significant  improvements. 


Character  of  Estuarine  Flows  and  Boundary  Stresses 

5.  Two  aspects  of  estuaries  require  a  method  of  evaluation  of  estuarine 
currents  different  from  the  usual  steady,  uniform  case.  First,  their  oscilla¬ 
tory  nature  creates  time -dependent  accelerations  that  alter  the  thickness  and 
velocity  profile  of  the  bottom  boundary  layer.  Second,  freshwater  runoff 
meeting  saline  ocean  water  creates  longitudinal,  and  usually  lateral  and  ver¬ 
tical,  density  gradients  that  cause  the  ambient  vertical  velocity  profile  to 
be  substantially  different  from  the  classic  forms.  Figure  1  illustrates  the 
general  nature  of  the  vertical  and  longitudinal  structures  of  estuarine  flows. 


( _ "  INDICATES  THAT  DIRECTION  OF 

CURRENT  REVERSES  WITH  TIOAL 
PHASE,  BUT  THAT  PRE'  OMINANT 
CURRENT  IS  IN  DIRECTION  OF 
LONGER  ARROW 

Figure  1.  Schematic  vertical  and  longitudinal  current 
structure  of  estuarine  flows 


6.  Unlike  unidirectional  flows,  in  which  the  entire  depth  of  flow  may 
be  occupied  by  the  boundary  layer,  estuarine  boundary  layers  grow  between  suc¬ 
cessive  slack- water  periods;  thus,  their  thickness  may  be  small  with  respect 
to  the  water  depth  duri.ig  a  significant  portion  of  the  tidal  cycle.  This  has 


significant  implications  with  regard  to  use  of  standard  velocity  profile  ecua- 
tions  and  constant  coefficients  for  computing  stresses.  Several  researchers 
have  observed  that  shear  stresses  are  dependent  on  the  pressure  gradient  as 
well  as  the  flow  speed.  Gordon  (1975)  measured  Reynolds  stresses  in  a  tidal 
current  and  concluded  that  bed  shear  stresses  were  considerably  larger  during 
decelerating  flows,  i.e.,  when  the  pressure  gradient  was  in  opposition  to  the 
direction  of  the  current,  than  during  accelerating  flows. 

Soft  Sediment  Boundaries 

7.  The  bed  elevation  in  open  channels  is  usually  defined  as  the  loca¬ 
tion  at  which  the  flow  velocity  is  zero.  Where  beds  are  rigid,  this  defini¬ 
tion  is  appropriate;  but  for  mobile  sand  beds  the  definition  lacks  precision 
and  for  soft  muds  it  may  be  substantially  erroneous. 

8.  In  many  estuarine  situations,  sediment  concentration  and  bulk  den¬ 
sity  increase  rapidly  near  the  bed,  and  often  do  so  without  a  clearly  defined 
interface  between  water  and  mud.  In  the  words  of  W.  R.  Parker  (Blackdown 
Consultants,  Taunton,  England),  as  the  bed  is  approached,  the  suspension 
changes  gradually  "from  muddy  water  to  watery  mud  to  bed"  without  a  definite 
interface.  Under  these  conditions,  defining  the  bed  elevation  becomes  a  mat¬ 
ter  of  selecting  an  arbitrary  definition,  usually  based  on  density. 

9.  A  density-defined  bed  interface  need  not  correspond  to  the  zero 
velocity  location,  since  at  any  given  moment  the  high-density  suspension  may 
flow  to  a  significant  depth  or  it  may  possess  sufficient  structure  (stationary 
suspension)  to  hold  it  apparently  motionless  In  either  case,  ordinary 
definitions  of  the  boundary  and  boundary  layer  are  often  inappropriate. 

Scone 

10.  This  work  is  intended  to  be  limited  to  conditions  characteristic  of 
placement  of  dredged  material  in  estuarine  sites  in  the  United  States,  though 
the  applicability  may  be  wider  than  that.  The  emphasis  is  on  sites  where  the 
transport  mechanisms  are  primarily  currents  with  some  contribution  by  non¬ 
breaking  waves.  A  parallel  investigation  on  wave -dominated  environments  is 
reported  separately  (Madsen  and  Wikramanayake  1988).  This  report  is  concerned 
with  existing  theory  and  experiments  and  how  they  may  be  best  adapted  to 
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engineering  calculations  for  sediment  transport  at  open-water  dredged  material 
placement  sites.  Sediment  transport  calculations  will  be  dealt  with 
separately. 
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PART  II:  HYDRODYNAMIC  AND  BOUNDARY  LAYER  FUNDAMENTALS 


Basic  Equations 

11.  Present  hydrodynamic  analyses  almost  universally  begin  with  the 
Eulerian  view  Navier-Stokes  equations  for  conservation  of  mass  and  momentum  in 
free  surface  flows,  which,  expressed  in  tensor  notation,  are 


Dp 

Tix" 


dp 

It 


0 


(1) 
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3xl 
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(2) 


where* 

p  -  fluid  density 
x  -  spatial  coordinate 
t  -  time 

u  -  flow  velocity 
P  -  pressure 

g  -  acceleration  due  to  gravity 
r  -  bed  shear  stress 

i,j,k  -  coordinate  axis  designator,  1,  2,  or  3  with  repeated  index 
summation 

-  indicates  instantaneous  value 

An  equation  for  conservation  of  constituents  in  transport  is  added 


P 


D9n 

TfiT 


d 


%)  +  K 


(3) 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  A).  The  notation  for  Parts  II  and  III  is  listed 
separately  from  that  for  Part  IV,  as  symbols  may  have  different  meanings  for 
the  different  areas  of  study  covered  in  each  Part. 


where 


6  -  concentration  of  constituent 

n  -  subscript  indicating  the  constituent  (e.g.  salinity, 
temperature,  or  sediment) 

E^  -  Fickian  diffusion  coefficient 

4>  =  sum  of  sinks  and  sources  for  8 


An  equation  of  state  completes  the  set  of  equations: 


p  -  /(©n> 


(4) 


For  turbulent  flows  in  large  bodies  of  water,  the  instantaneous  variables  of 
Equations  1-4  are  replaced  with  a  summation  of  a  time-mean  value  (averaged  on 
the  order  of  minutes)  and  an  instantaneous  fluctuation  about  that  mean;  for 
example , 

-  (  S  'S 

u  -  u  +  u'  1  ' 


where  the  prime  indicates  a  quantity  fluctuating  in  time.  Replacement  of  the 
instantaneous  values  with  the  form  of  Equation  5  and  performing  algebraic 
manipulation  results  in  number  of  additional  terras.  Performing  a  second 
time -averaging  step  v which  causes  some  terms  to  become  equal  to  zero), 
neglecting  some  terms  that  make  apparently  minor  contributions  (including 
spatial  variations  in  density  except  in  the  pressure  term),  and  adding  a  term 
for  the  Coriolis  effect  results  in 
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p  -  /(e„) 
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where 


e  -  alternating  tensor 
0  -  Coriolis  effect  parameter 

n  -  index  for  multiple  constituents  in  transport 

12.  The  overbar  terms  in  equation  set  6  are  time  averages  of  the  prod¬ 
ucts  of  fluctuating  variables  and  are  described  as  Reynolds  stresses  (Equa¬ 
tion  6b)  or  Reynolds  fluxes  (Equation  6c) .  Defining  these  terms  becomes  the 
turbulence  closure,  which  is  a  fundamental  problem  in  practical  use  of  the 
equations.  That  problem  is  discussed  in  Part  III  of  this  report. 

Laminar  and  Turbulent  Flows 

13.  Flows  of  interest  here  are  usually  fully  turbulent  in  terms  of  the 
overall  Reynolds  number  criterion.  Boundary  layers  are  not  necessarily  fully 
turbulent,  though.  The  boundary  Reynolds  number  of  Nikuradse  has  been  adopted 
by  marine  researchers,  including  Sternberg  (1968),  who  defined  the  turbulence 
characteristics  as  follows: 


a . 

Smooth  turbulent  when 

Rp  <  5.5 

b. 

Transitional  when  5.5 

<  Rp  <  165 

c . 

Rough  turbulent  when 

165  <  Rp 

where 


Rp 


u.ks 


(7) 


and 

u*  -  shear  velocity 
ks  =  roughness  element  size 
v  -  kinematic  viscosity 


Boundary  Layer  Structure 

14.  Boundary  layers  in  open  channel  flow  are  the  zones  near  the  boundary 
where  drag  exerted  by  the  bed  has  a  pronounced  effect  on  the  flow.  The 
boundary  layer  is  characterized  by  steep  velocity  gradients  normal  to  the  bed 
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and  a  corresponding  increase  in  turbulent  and  viscous  stresses. 

15.  The  boundary  layer  thickness  can  be  defined  in  several  ways.  In 
steady  viscous  flow  it  is  commonly  defined  as  the  point  at  which  the  velocity 
reaches  99  percent  of  the  free  stream  velocity  or  a  similar  measure  in  terms 
of  the  inviscid  flow  condition.  While  none  of  these  definitions  is  par¬ 
ticularly  precise,  they  are  useful  in  describing  the  flow  profiles  and  bed 
stresses . 

16.  Turbulent  boundary  layers  exhibit  zones  in  which  properties  of  the 
flow  differ.  Figure  2  illustrates  zone  subdivisions.  Figure  2a  shows  a  prev¬ 
alent  subdivision  of  steady- flow  boundary  layers.  The  size  and  existence  of 
the  viscous  (laminar)  sublayer  and  buffer  zone  are  affected  by  the  degree  of 
turbulence  and  texture  of  the  boundary,  and,  for  fully  rough  turbulent  flow, 
they  may  disappear. 

17.  For  smooth  turbulent  flows,  the  viscous  sublayer  thickness  is  much 
greater  than  ks  ,  while  for  transitional  and  rough  turbulent  flows  it  is 
approximately  equal  to  and  much  less  than  k,  ,  respectively.  Thus  as  tur¬ 
bulence  increases,  bottom  roughness  plays  an  increasing  role  in  the  nature  of 
the  boundary  layer. 

18.  Figure  2b  shows  the  structure  terminology  used  by  Soulsby  (1983) 

for  marine  boundary  layers.  The  bed  layer  can  be  either  laminar  (when  the  bed 
is  sufficiently  smooth)  or  turbulent.  Soulsby  refers  to  the  second  zone  as 
the  log  layer,  where  classical  profile  equations  are  thought  to  apply,  and  to 
the  third  zone  as  the  outer  layer.  In  the  outer  layer,  structure  depends  on 
the  free  streamflow  behavior,  and  thus  the  classical  equations  do  not  apply. 
Above  the  outer  layer  is  the  free  stream  zone,  where  velocity  profile  is 
essentially  unaffected  by  the  bed. 

19.  The  boundary  layer  structure  in  estuarine  flows  is  in  general  sig¬ 
nificantly  altered  by  the  presence  of  both  tidal  currents  and  wind- generated 
surface  waves.  The  superposed  current  and  high-frequency  oscillatory  water 
particle  motion  beneath  a  water  wave  interact  in  a  nonlinear  manner  and  result 
in  a  change  in  both  the  magnitude  and  direction  of  the  bed  stresses.  For  com¬ 
bined  wave  and  current  flows,  a  wave  boundary  layer  develops  in  the  zone  imme¬ 
diately  above  the  bottom.  The  thickness  of  this  boundary  layer  is  given  by 
xu.ew/(7  ,  where  k  is  von  Karman’s  turbulence  constant,  a  is  wave  fre¬ 
quency,  and  u„CH  is  wave-  and  current- induced  shear  velocity.  Within  the 
wave  boundary  layer,  turbulence  production  is  due  to  both  wave  and  current 


b.  Soulsby  (1983)  subdivision  of  marine  boundary  layers 
Figure  2.  Boundary  layer  subdivisions 


motion.  A  large-scale  rotating  current  boundary  layer  extends  from  the  bottom 
to  the  free  stream  zone  (when  the  latter  is  present)  or  the  water  surface.  In 
the  portion  of  the  current  boundary  layer  that  extends  above  the  wave  boundary 
layer,  turbulence  is  associated  with  the  oscillatory  tidal  current  only. 

Boundary  Stresses 

20.  Stresses  experienced  by  the  sediment  bed  are  complex  products  of  the 
interaction  of  flow  anti  bed.  For  convenience,  these  stresses  are  commonly 
divided  into  tangential  and  normal  stresses.  Many  analyses  divide  them  into 
surface  and  form  stresses.  Surface  stress  (or  drag)  is  that  which  would  occur 
if  the  bed  were  plane,  and  it  can  be  described  in  terms  of  the  height  by  which 
sediment  grains  protrude  into  the  flow.  Form  stress  is  that  caused  by 
irregularities  in  the  boundary,  such  as  ripples  or  dunes. 

21.  Bed  stresses  contributing  to  sediment  transport  are  often  conceived 
of  and  described  as  tangential,  even  when  they  contain  normal  components  (as 
in  form  drag).  In  some  analyses,  the  normal  lift  force  is  assumed  to  be  pro¬ 
portional  to  the  tangential  shear  stress  and  thus  is  included  implicitly. 

These  approaches  cause  difficulties  in  estuarine  and  coastal  environments, 
where  soft  mud  beds  can  deform  under  normal  stresses,  changing  the  relation¬ 
ship  between  tangential  and  normal  stresses,  or  where  porous  beds  can  experi¬ 
ence  large  (and  fluctuating)  normal  stresses  and  flow  through  the  bed. 

22.  For  these  reasons  and  the  soft  beds  problems  cited  in  para¬ 
graphs  7-9,  shear  stress  exerted  on  the  sediment  bed  by  the  flow  is  an 
extremely  unsatisfactory  parameter  for  sediment  transport  calculations. 
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PART  III:  BOUNDARY  LAYERS  AND  STRESSES  DUE  TO  CURRENTS 


23.  Use  of  equation  set  6  involves  selecting  a  turbulence  closure  to 
represent  the  Reynolds  stresses  and  fluxes  and  making  simplifying  assumptions 
to  the  degree  necessary  for  the  solution  method  used  and  suitable  for  the  pur¬ 
pose.  Computational  difficulty  and,  more  recently,  computer  cost  have  forced 
use  of  rather  extreme  simplifications  to  obtain  estimates  of  shear  stresses 
and  velocities,  in  spite  of  the  loss  of  realism. 

24.  To  gain  an  appreciation  for  how  simplification  of  the  equations 
affects  the  end  result,  the  following  results  are  presented. 

25.  Equation  6b  can  be  integrated  over  width  so  that  spatial  variations 
are  considered  only  in  the  vertical  and  longitudinal  directions.  Integrating, 
dropping  the  Coriolis  effect,  and  splitting  the  pressure  term  into  a  surface 
slope  term  and  a  density  gradient  term  yields  a  two-dimensional  (vertical) 
equation  of  motion 


ABC  D 


where 

u  -  time -mean  velocity 

x  -  coordinate  direct5.on  variable 

r)  -  displacement  of  water  surface  from  mean  water  level 
tsw  -  sidewall  shear  stress 

'  -  an  instantaneous  turbulent  fluctuation  of  a  variable 
and  A  is  the  unsteadiness  term  (temporal  acceleration),  B  is  the 
nonuniformity  terra  (convective  acceleration),  C  is  the  surface  slope  term, 
D  is  the  longitudinal  density  gradient  term,  and  E  and  F  are  shear 
stresses  arising  from  boundary  resistance. 
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26.  Evaluation  of  the  relative  importance  of  the  terms  in  Equation  8  can 
be  accomplished  by  an  order  of  magnitude  analysis  or  by  use  of  experimental 
data.  Abbott  (1960)  did  the  former  by  assuming  that  the  tide  could  be 
described  as  a  simple  progressive  wave.  For  parameters  typical  of  the  Thames 
River  he  found  that  the  unsteady  and  surface  slope  terms  (A  and  C)  were  the 
most  important,  followed  by  nonuniform  flow  and  density  gradient  terms  (B 

and  D) . 

27.  An  experimental  evaluation  of  the  terms  in  Equation  8  was  performed 
using  data  from  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  salin¬ 
ity  flume  (test  14)  in  which  a  10-min  tide  of  0.10-ft*  amplitude  was  induced 
in  the  ocean  headbay  (salinity  30  ppt)  while  a  freshwater  flow  of  0.0075  cfs 
was  introduced  at  the  upstream  end  of  the  327-ft-long,  0.75-ft-wide  flume. 

The  vertical  salinity  gradient  in  this  test  was  moderate,  with  about  a  30  per¬ 
cent  difference  between  salinities  at  the  bottom  and  surface  at  a  point  40  ft 
from  the  ocean  end  of  the  flume . 

28.  Using  a  three-point  piecewise  linear  fit  (e.g.,  see  Halliwell  and 
O'Connor  1968)  to  data  points  in  the  vertical  and  difference  forms  for  the 
derivatives,  terms  A-D  in  Equation  8  were  evaluated  at  one  station  in  the 
salinity  flume.  Figure  3  shows  the  individual  terms'  variation  in  magnitude 
plus  their  sum  over  a  complete  cycle.  At  the  bottom  of  the  figure  the  depth- 
averaged  current  velocity  and  water-surface  elevation  histories  at  that 
station  are  shown.  Terms  A-D  were  nondimens ionalized  by  dividing  by  pgh  , 
where  h  is  the  mean  water  depth. 

29.  As  in  Abbott's  analysis,  the  surface  slope  and  unsteady  terms  were 
the  largest,  but  at  times  during  the  cycle  the  density  and  nonuniform  flow 
terms  became  significant  in  comparison.  The  nonuniform  flow  term  should  not 
be  expected  to  be  large  since  the  flume  is  prismatic.  Using  the  three -point 
piecewise  fit  and  considering  variation  of  p  with  respect  to  y  (vertical 
coordinate)  resulted  in  about  a  5  percent  reduction  in  the  magnitude  of  the 
surface  slope  terms  for  the  flume  data.  This  suggests  that  for  similar  condi¬ 
tions  that  can  be  adequately  described  by  linear  approximations,  the  variation 
of  p  with  y  may  be  neglected  for  individual  terms  and  the  total  shear 
except  when  the  direction  of  flow  is  reversing.  However,  these  analyses 


*  A  table  of  factors  for  converting  non-Sl  units  of  measurement  to  SI 
(metric)  units  is  found  on  page  4, 
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showed  that  it  may  not  be  appropriate  to  automatically  assume  further 
simplifications  to  Equation  8  for  all  cases.  Rather,  for  each  application  an 
evaluation  must  be  made  as  to  the  proper  simplifying  assumptions. 

30.  Figure  4  shows  how  the  actual  slope  is  approximated  by  a  typical 
steady- flow  bed  shear  stress  equation  of  the  form 


t  - 


(9) 


where  fc  is  a  constant  friction  factor  with  an  empirically  derived  steady- 
flow  value  of  0.02  for  the  salinity  flume,  and  ux  is  a  flow  velocity  at  a 
specified  distance  above  the  bed. 

31.  Figures  3  and  4  show  that,  for  the  fairly  simple  case  of  stratified 
flow  in  a  prismatic  flume,  using  the  surface  slope  or  a  simple  velocity 
squared  equation  for  calculating  the  boundary  stresses  can  lead  to  signifi¬ 
cantly  low  magnitude  estimates.  At  small  values,  even  the  direction  of  stress 
can  be  in  error.  In  more  complicated  flows  characteristic  of  natural  environ¬ 
ments,  the  errors  can  be  even  larger. 

Steady.  Uniform  Flows 

32.  For  steady,  constant -density,  uniform  flows,  the  coefficient  fc  in 
Equation  9  is  equal  to  one  fourth  the  Darcy  pipe  flow  friction  factor  f  ,  and 
when  the  flow  is  also  rough  and  turbulent,  it  can  be  related  to  the  Chezy  C  , 
Manning's  n  (non-SI  units),  and  effective  roughness  ks  coefficients  by 


fc  - 


2g 

C2 


2gn2 

(l .  49R1/6)2 


(10) 


where  R  is  the  hydraulic  radius. 

33.  Shear  stress  can  be  calculated  by  Equation  9  with  an  estimated  value 
of  fc  ,  or  by  fitting  a  theoretical  profile  to  observed  velocity  variation  in 
the  vertical.  Table  1  shows  the  commonly  used  equations  in  a  form  that  pro¬ 
vides  an  estimate  of  the  roughness  coefficient.  The  Prandtl  equations  are 


19 


0.6 


Figure  4.  Flume  boundary  stress  by  equation  of  motion  and 
by  constant  friction  factor 


CURRENT  VELOCITY,  FPS 


limited  to  the  near-wall  portion  of  the  boundary  layer  because  of  Prandtl's 
assumption  that  shear  stress  in  the  fluid  is  not  a  function  of  distance  from 
the  wall.  That  restriction  is  often  ignored,  but  using  velocities  more  than 
15  percent  of  6  (the  boundary  layer  thickness)  away  from  the  wall  can  lead 
to  shear  stresses  that  are  too  high  by  significant  amounts  (Lyles  and  Woodruff 
1972).  The  corresponding  velocity  profile  equations  are  as  follows: 

a.  Modified  von  Karman- Prandtl 


u 


-1  In 
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+  0.0338 
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(11) 


b.  Velocity  defect 
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y  >  0.156 

2.5 

y  <  0.155 

0 

y  >  0.155 

c .  Power  law 


JL  =8.74 
u. 


(13) 


where 

y  -  distance  above  the  bed 
U0  -  free-stream  steady  velocity 
6  -  boundary  layer  thickness 
Cyi ,  Cyj  *  constants 
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Unsteady.  Nonuniform  Flows 

34.  The  challenge  in  solving  for  turbulent,  unsteady,  nonuniform  flows 
lies  in  finding  a  suitable  way  to  solve  for  the  Reynolds  stress  and  flux  terms 
of  Equation  6.  In  the  simplest  form,  the  terms  are  integrated  over  space  (as 
in  Equation  8)  and  an  expression  such  as  Equation  9  for  boundary  shear 
stresses  is  employed.  A  somewhat  more  complicated  eddy  viscos _ty  approach  is 
used  for  internal  flow  stresses  and  sometimes  for  boundary  stresses  as  well. 
(An  analogous  approach  is  used  for  fluxes.)  Progressively  more  complex  (and 
sometimes  more  realistic)  descriptions  of  the  Reynolds  terms  are  attempted 
through  turbulence  transport  models,  sometimes  referred  to  as  higher  order 
turbulence  closures. 

35.  The  array  of  possible  solutions  to  the  Reynolds  terms  spans  a  con¬ 
tinuum  of  complexity  and  physical  realism.  Here,  that  continuum  is  broken  for 
convenience  into  steady  uniform  flow  equations,  eddy  viscosity  models,  and 
turbulence  transport  models.  The  laminar  case  is  first  briefly  described, 
then  turbulent  flows  are  considered  in  more  detail. 

Laminar  flow 

36.  Schlicting  (1968)  solved  an  approximate  expression  for  laminar 
boundary  layer  flows  that  occur  when  an  oscillating  current  occurs  adjacent  to 
a  no-slip  boundary.  Although  estuarine  flows  are  rarely  laminar,  they  may  be 
so  very  near  the  boundary  and  near  slack  water.  This  simplified  case  illus¬ 
trates  the  behavior  of  oscillatory  boundary  layers  and  gives  a  starting  point 
for  further  development. 

37.  In  Schlicting' s  solution  a  free-stream  velocity  of 


u  -  UD  sin  o 


(14) 


where 

UQ  =  maximum  free-stream  velocity 
a  -  2tt/T,  frequency  of  oscillation 
T  «*  period  of  oscillation 

adjacent  to  a  fixed,  no-slip  boundary  results  in  a  near-boundary  velocity  pro¬ 
file  given  by 
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u(y,t)  -  UQ  sin  at  [l  -  e  (cos  0y  sin  at  -  sin  j3y  cot  at  )] 


<1S) 


where 


0  •  jo /2a 


Equation  15  is  a  damped  harmonic  function  in  y  ,  the  distance  away  from  the 
boundary-  Figure  5  illustrates  the  form  of  the  profile  at  various  times  dur¬ 
ing  T  .  Note  that  since  the  boundary  layer  velocity  overshoots  the  free- 
stream  velocity,  a  standard  definition  of  steady-flow  boundary  layer  thick¬ 
ness,  such  as  Che  distance  between  the  wall  and  the  point,  at  which  u  -  0.99Uo 
does  not  have  the  same  physical  significance  in  the  oscillatory  case.  If  the 
boundary  layer  is  defined  as  that  zone  where  the  velocity  differs  signifi¬ 
cantly  from  the  free-stream  velocity,  the  point  at  which  0y  -  *  is  a  more 
appropriate  measure.  Thus,  the  oscillatory  boundary  layer  thickness  would  be 
expressed  as 


6  -*J  »  jjvT  (16) 

For  laminar  flows  of  tidal  periods,  this  yields  a  maximum  thickness  of  1.5  ft. 
Turbulent  flows 

38.  The  solution  for  laminar  flow  (Equation  15)  does  not  extend  to 
boundary  layers  in  turbulent  flows.  Many  researchers  have  assumed  or  at¬ 
tempted  to  show  that  the  boundary  layer  in  oscillatory  flow  of  tidal  periods 
constitutes  a  considerable  portion  of  the  depth  of  flow.  Sternberg  (1968) 
measured  Puget  Sound  tidal  velocities  at  three  points  within  1  m  of  the  bed 
and  found  that  velocity  profiles  at  sites  with  intense  flow,  complex  geometry, 
and  small  density  gradients  could  be  satisfactorily  represented  (the  curve  fit 
the  data  points  with  some  scatter)  by  a  logarithmic  profile  62  percent  of  the 
time.  Since  the  logarithmic  portion  of  the  velocity  profile  (see  para¬ 
graph  33)  is  limited  to  about  the  lower  15  percent  of  the  boundary  layer,  a 
1-m- thick  log  layer  implies  a  boundary  layer  about  5-10  m  thick. 
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Lesht  (1978),  measuring  the  frequency  of  turbulent  fluctuations  and  mean  flow 
speed  in  the  New  York  Bight,  found  that  velocities  within  1  m  of  the  bed 
demonstrated  varying  degrees  of  correlation  with  a  log  profile,  depending 
largely  on  the  time  period  over  which  velocities  were  averaged.  For  averaging 
periods  greater  than  4  min,  the  correlation  coefficients  ranged  from  0.55  to 
0.90.  Lesht  attributed  low  correlation  at  shorter  averaging  periods  to 
wind/wave-  induced  velocity  fluctuations. 

39.  Mehta  and  Christensen  (1977)  applied  the  modified  von  Karman- Prandtl 
equation  (Equation  11)  to  flow  in  tidal  inlets  using  several  velocity  measure¬ 
ments  within  1  m  cf  the  bed.  They  reported  a  satisfactory  fit  to  the  data 
points  by  the  theoretical  velocity  profile.  Ludwick  (1973),  assuming  the 
boundary  layer  thickness  to  be  equal  to  the  flow  depth,  applied  the  velocity 
defect  concept  (Equation  12)  to  Chesapeake  Bay  flows  with  reasonable  results. 

40.  The  presence  of  vertical  or  horizontal  density  gradients  further 
complicates  consideration  of  estuarine  velocity  profiles.  Density  gradients 
distort  the  velocity  profile  from  that  of  well-mixed  flows,  resulting  in  a 
tidally  averaged  circulation  pattern  like  that  of  Figure  1.  At  times  the 
near-bed  velocity  may  be  in  the  opposite  direction  from  that  of  the  free 
stream,  due  either  to  phase  shifts  as  shown  in  the  laminar  case  or  to  the 
pressure  gradient  caused  by  density  differences.  Even  if  we  may  marginally 
assume  1-m-thick  logarithmic  profiles  in  well-mixed  tidal  flows,  such  an 
assumption  is  quite  often  invalid  for  partially  or  strongly  stratified  flows 
as  shown  in  Figure  1 . 

Using  steady,  uniform  flow  equations 

41.  Sternberg  (1968),  using  a  log  profile,  calculated  values  of  fc 
(Equation  9)  for  a  characteristic  velocity  1  m  above  the  bed.  He  found  a  mean 
value  of  fc  -  0.0062  ,  with  95  percent  confidence  limits  of  0.0017  and  0.022. 
For  Reynolds  numbers  greater  than  1.5  x  105,  results  were  less  variable  and 
the  mean  value  of  fr  was  0.0060. 

42.  Bowden,  Fairbairn,  and  Hughes  (1959)  found  that  reasonable  shear 
stresses  resulted  from  Equation  9  using  a  value  of  fc  of  0.007  and  a  charac¬ 
teristic  velocity  measured  at  a  distance  above  the  bed  equal  to  5  percent  of 
the  water  depth. 

43.  Though  several  of  these  researchers  have  found  that  velocity  pro¬ 
files  can  often  be  well  approximated  by  a  logarithmic  profile,  it  does  not 
necessarily  follow  that  the  assumption  of  a  log  profile  yields  a  satisfactory 
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estimate  of  shear  stress.  There  is  usually  a  great  deal  of  noise  in  velocity 
measurements ,  and  any  continuous  (in  space)  expression  is  generated  on  a  best- 
fit  basis  to  that  data.  Thus,  a  theoretical  velocity  profile  may  "look"  like 
a  satisfactory  approximation  to  the  measured  velocity  points  without  correctly 
predicting  the  shear  stress  at  the  bed. 

44.  Table  2  lists  values  of  fc  obtained  by  several  investigators  and 
those  resulting  from  Equations  10,  11,  and  12.  Although  there  is  a  consider¬ 
able  range  between  the  highest  and  lowest  values,  most  are  equal  to  about 
0.007,  As  previously  mentioned,  part  of  the  difficulty  in  defining  fc  from 
field  data  is  that  it  is  time  dependent,  in  that  flow  may  be  changing  from 
hydraulically  smooth  to  hydraulically  rough,  the  bed  itself  may  be  changing 
with  time  and  space,  and  the  pressure  gradient  dependency  of  friction  (Ludwick 
1973;  Smith  and  McLean  1977)  is  not  included.  Under  these  conditions,  average 
shear  stresses  may  be  computed  from  suitably  time-  and  space -averaged  velocity 
profiles,  but  instantaneous  local  shear  stresses  may  differ  substantially  from 
such  averages.  It  is  probable  that  the  field-determined  values  of  fc  are 
more  appropriate  for  total  shear  than  for  bed  shear,  particularly  for  noncohe- 
sive  bed  cases.  The  values  in  Table  2  suggest  that  fc  for  the  flat  bed  case 
may  be  as  much  as  an  order  of  magnitude  less  than  field-derived  values.  In 
the  absence  of  better  information,  one  of  the  forms  of  Equation  10  can  be  used 
for  a  constant  friction  coefficient,  and  an  estimate  of  a  time-varying  bed 
surface  shear  stress  can  be  obtained  by  multiplying  the  values  from  Equa¬ 
tion  10  by  the  cube  root  of  the  ratio  of  kg  based  on  grain  size  to  an 
effective  total  k5  . 

Eddy  viscosity  models 

45.  Next  in  complexity  after  the  steady- flow  equations  come  the  simple 

turbulence  closure  models.  They  are  typically  based  on  a  turbulent  eddy  vis¬ 
cosity  coefficient  that  is  analogous  to  the  kinematic  viscosity  v  and 

is  used  in  an  expression  for  shear  stress 


r 


(17) 


which  parallels  the  expression  for  laminar  shear  stress.  Knight  (1978)  and 
Sleath  (1984)  provide  excellent  reviews  of  these  models. 
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46.  The  difficulty  with  this  approach  is  defining  an  appropriate  value 
for  the  eddy  viscosity.  Since  it  is  a  property  of  the  flow  instead  of  the 
fluid  (like  kinematic  viscosity),  eddy  viscosity  is  not  constant  in  space  or 
time.  In  the  simplest  case  it  is  assumed  constant.  In  progressively  more 
realistic  (and  computationally  more  complex)  approaches  it  is  assumed  to  be 
constant  over  two  or  three  horizontal  layers  in  the  water  column.  Only  one 
such  model  is  presented  here. 

47.  Three- layer  model.  Kaj iura  (1968)  proposed  a  widely  cited  three- 
la_yei:  model  that  defines  an  inner  layer  next  to  the  bed,  an  overlap  layer,  and 
an  outer  layer.  The  eddy  viscosity  for  each  layer  is  given  by 

(u  for  smooth  flow  (18) 

‘'t.  inner  “  utl  *  1  2 . 71  KU ,6  y 


•'t. overlap  "  m  *u.y 


(19) 


^t. outer  "  ut3  m 


(20) 


where 

k  -  0.4 

-  inner  layer  thickness 
60  -  thickness  of  all  three  layers 

and 


*o-0.05^ 


f  12v 

u. 


smooth  beds 


rough  beds 


The  overlap  layer  disappears  if  >  Sa  . 

48.  For  an  oscillatory  free-stream  velocity  of  the  form 


(21) 


(22) 
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(23) 


u  »  UQ  cos  at 


the  equation  of  motion  yields  for  the  boundary  layer  velocity 


u  =  REAL  UD  (l  -  e  lY  ♦  sinh  Ojy)e1<rt 


(24) 


where 


Ai 


(<*x  ♦  a2)e°l4° 

cosh  (ai^i)  ♦  a2  sinh  (a1i1) 


(25) 
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—  for  smooth  beds 
v 


To 

0.185  icu.^ 


for  rough  beds 
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(26) 


(27) 


and 

i=/T 


From  these  equations  it  follows  that  the  bed  shear  stress  is 


r 


=  i/tREAL  [axU^Ai  cosh  axy  - 


(28) 


Computation  of  velocity  profiles  and  stresses  using  Equations  18-28  requires 
knowledge  of  the  free-stream  velocity,  oscillation  period,  and  effective 
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roughness  size.  Since  the  equations  also  include  the  shear  velocity,  an 
iterative  solution  is  required. 

49.  Mixing  length  model.  Mixing  length  models  return  to  the  Prandtl 
hypothesis  that 


(29) 


and  rely  on  one  of  several  expressions  for  the  value  of  i  ,  the  mixing 
length.  One  set  of  expressions  commonly  used  is 


40t/ 

0  <  y  <  yx 

xy  - 

u. 

a,  6 

yx<y<  6 

(30) 


where  a,  and  yi  are  empirical  parameters. 

50.  Continuous  analytic  forms.  More  recently,  eddy  viscosity  formula¬ 
tions  in  width- integrated  and  three-dimensional  numerical  models  have  been 
expressed  as  a  continuous  function  over  depth  of  a  constant  form  (e.g., 
parabolic)  with  a  magnitude  expressed  as  a  spatially  or  temporally  varying 
coefficient,  such  as  the  following  (King  1988): 


utxy  ~  vtR 


( 


B  + 


y 

TT 


■<  “ 

J  - 


(31) 


where 

i/txy  «  spatially  and  temporally  varying  eddy  viscosity 
utR  -  reference  value  of  eddy  viscosity 
A,B,C  -  constants 

and  i/yj  is  modified  by  an  appropriate  factor  in  stably  stratified  flows. 
51.  Reasonably  good  agreement  with  field  data  has  been  obtained  in 
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San  Francisco  Bay  using  this  general  approach  (Cheng,  Wang,  and  Gartner  1988). 
Equation  31  has  been  used  by  King  (1988)  in  the  three-dimensional  numerical 
model  RMA-10,  which  has  been  applied  to  estuarine  flow  problems  (Richards  and 
Bach  1988;  McAnally,  Letter,  and  Thomas  1986)  with  satisfactory  results. 
Advantages  and  disadvantages 

52.  These  steady- flow  and  simple  eddy  viscosity  approaches  offer  the 
advantages  of  simplicity  and  relatively  small  computational  effort.  They  are 
also  entirely  adequate  for  many  estuarine  applications,  including  water  level 
and  horizontal  circulation  calculations  and  even  depth- integrated  sediment 
transport  where  deposition  rates  are  of  paramount  interest.  They  are  less 
useful  and  sometimes  possibly  misleading  when  used  for  vertical  circulation, 
near-bed  velocity  profiles,  and  sediment  erosion  problems.  These  latter  cases 
require  that  the  simpler  approaches  be  used  with  extreme  caution  or  that  more 
general  methods  be  used. 

53.  The  eddy  viscosity  approaches  can  be  modified  for  some  effects  (such 
as  stratification)  by  analytical  and  empirical  methods.  These  methods  are 
described  later  in  this  part. 

Stratification  Effects 

54.  Methods  of  correcting  for  density  stratification  (e.g.,  due  to  tem¬ 
perature,  salinity,  and/or  suspended  sediment)  in  determination  of  boundary 
shear  and  velocity  profiles  in  surface  water  flows  are  necessary  if  the 
turbulence  calculation  method  does  not  explicitly  account  for  it. 

55.  French  (1977)  found  that  in  stratified  flow  the  Chezy  coefficient  C 
is  a  function  of  the  Reynolds  number,  boundary  roughness,  and  density  gradi¬ 
ent.  He  started  with  the  basic  relationship  between  C  and  the  shear 
velocity,  u*  , 


C  = 


(32) 


where  u  is  the  average  flow  velocity.  Then  he  derived  the  following  equa¬ 
tion  for  the  vertical  velocity  gradient  in  stratified  flows  using  Prandtl's 
mixing  length  assumption: 
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(33) 


du  _  pu? 

-  2  2  2  dp 

icy  pu.  -  afic2gyzig 

k  ^  j 


where  French  takes  a£  -  5.0  .  Integration  of  this  equation  gives  the  verti¬ 
cal  velocity  profile  u(y)  as  a  function  of  the  vertical  density  profile 
p(y)  .  Incorporating  Equation  33  into  Equation  32  gives 

C4|^)dy  (34) 

y<3 


where  yQ  is  the  bed  roughness  height.  Accounting  for  the  effect  of  strati¬ 
fication  on  the  velocity  profile  requires  solving  Equations  32,  33,  and  34  for 
a  given  density  profile  for  u*  and  iteratively  for  C  and  u(y)  .  Compari¬ 
son  of  C  predicted  using  Equation  34  with  a  few  experimental  data  sets 
indicated  good  agreement. 

56.  McCutcheon  (1981),  following  an  approach  similar  to  that  of  French, 
started  with  the  von  Karman-Prandtl  law  of  the  wall,  incorporated  the  similar¬ 
ity  analysis  of  Monin  and  Obukhov,  and  developed  the  following  formula  for  the 
vertical  velocity  distribution  for  unidirectional,  hydraulically  rough,  con¬ 
tinuously  stratified  open  channel  flow: 


u(y)  .  1  | 

r  dy 

u«  *  J 

1 

°  y 

1  ♦  agKZy2 

’ - ^  ■ 

/( y) 

- 

(35) 


where 

/(y)  -  dp/d y 

a  -  ra£  with  r  the  ratio  of  mass  mixing  length  to  momentum  mixing 
length 

pm  -  depth- averaged  density 

Equation  35  can  be  substituted  into  the  integrand  of  Equation  34,  and  then 
solved  for  C  and  u*  in  the  same  manner.  McCutcheon  performed  several 
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laboratory  tests  in  which  salt  was  the  stratifying  agent  in  order  to  verify 
Equation  35.  One  limitation  of  the  methods  by  French  and  McCutcheon  is  that 
they  both  assume  the  law  of  the  wall  to  be  valid  over  the  entire  depth  of 
flow.  In  addition,  if  the  vertical  density  profile  is  not  known  a  priori, 
these  equations  for  the  velocity  profile  would  have  to  be  solved  coupled  with 
the  conservation  of  mass  equation  for  salt  or  sediment  (Equation  6c). 


Turbulence  Transport  Models 


57.  Recently,  unsteady- flow  research  has  shifted  away  from  use  of 
analytic  expressions  and  constants  for  eddy  viscosity/mixing  length  and  con¬ 
centrated  on  various  higher  level  approaches  that  consider  turbulence  charac¬ 
teristics  in  more  detail.  These  approaches  aim  to  develop  additional  expres¬ 
sions  based  on  constants  that  are  truly  constant  among  a  given  class  of 
problems  (e.g.,  open  channel  flow)  and  not  flow  dependent.  Two  basic  ap¬ 
proaches  are  given  in  the  literature  under  a  variety  of  titles.  Reviews  of 
the  several  turbulence  transport  approaches  are  available  (e.g.,  Markatos 
1986;  Task  Committee  on  Turbulence  Models  in  Hydraulic  Computations  1988); 
however,  only  two  are  presented  here. 

58.  Turbulence  transport  models  begin  with  the  Eulerian  conservation 
of  mass  and  momentum  equations  and  an  equation  of  state,  which  are  given  in 
Part  II  as  Equation  6. 

59.  In  turbulence  transport  models,  the  Reynolds  stress  and  flux  terms 
are  expressed  in  further  equations  that  involve  production,  transport,  and 
dissipation  of  turbulent  kinetic  energy.  Two  prominent  approaches  are  the 

k  -  e  model  and  the  level  2-1/2  model,  both  described  in  the  following 
paragraphs . 
k  -  e  model 

60.  In  k  -  e  models  the  turbulent  kinetic  energy  k  is  transported, 
produced,  and  dissipated.  These  models  have  been  widely  described  and  applied 
(e.g.,  see  Rodi  1982  and  1987).  The  conservation  of  turbulent  energy  equation 
is  expressed  as 

i  ^  *£]  * p-  * 0  -  ■  <36a> 
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where  Ps  represents  production  of  turbulence  by 

r 

p  _  „  3ui  3uJ  3ui 

8  "  I*;  1*1  -3xJ 


(36b) 


G  represents  buoyancy  contribution  to  turbulence, 


dp 

at 


(36c) 


where 

p0  -  reference  density 
<rt  -  turbulent  Prandtl  number 
€  -  dissipation  of  turbulence  expressed  by 

*«--««•£]  (36d) 

with  vt  equal  to  C„  (k2/0  •  The  notation  is  the  same  as  in  Part  II  and 
the  C's  are  constants.  If  the  equations  are  truly  general,  the  five  C's 
will  actually  be  constant  for  all  flows.  Their  values  are  given  by  Bradshaw, 
Cebeci,  and  Whitelaw  (1981)  and  Rodi  (1987)  as 

C„  -  0.09 

Cjjj  -  1.0 

Cel  -  0.77 
Ce2  -  1.45 
Ce3  -  1.9 

The  value  for  CeG  is  1  for  positive  G  and  0-0.2  for  negative  G  (Task 
Committee  1988). 

61.  The  k  -  e  model  offers  coefficients  that  are  more  reliably  con¬ 
stant  over  a  given  class  of  flows  than  are  the  empirical  coefficients  of 
simpler  methods  by  incorporating  more  physical  realism  at  the  cost  of  sharply 
increased  computational  effort.  While  the  transport  of  turbulent  energy  seems 
clearly  a  more  physically  reasonable  concept  than  simpler  eddy  viscosity 
models,  the  concept  of  transporting  dissipation  is  not  so  clearly  reasonable. 
Were  it  not  for  the  k  -  €  model’s  success  in  reproducing  a  number  of  complex 
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flow  fields  (Rodi  1987)  ,  the  concept  would  seem  suspect  because  of  the 
dissipation  equation. 

62.  More  serious  deficiencies  in  the  k  -  e  model  are  that  (a)  in  its 
basic  form  it  does  not  account  for  flow  stratification  effects  on  turbulence 
and  (b)  it  does  not  predict  the  flow  field  close  to  solid  boundaries  well 
without  special  adaptation  to  a  law  of  the  wall  form. 

63.  Takemitsu  (1988)  has  proposed  a  modification  to  the  k  -  £  model 
that  replaces  Equation  36d  with 
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which  increases  the  number  of  constants  to  nine.  He  specifies  the  range  of 
constant  values  as  shown  in  the  following  tabulation. 


Constant 


Minimum 

Value 


Probable 

Value 


Maximum 

Value 


Cki 

1 

3 

0(10) 

0.4Ckl 

0 • 5Ckl 

0.6Ckl 

Cu 

1 

1 

°(10) 

CfZ 

2 

Ce3 

2 

Cf* 

0 

1 

®(10) 

Ce5 

1 

c«  l 

°(10) 

Cf6 

-10 

0 

2 

C<7 

-5 

0 

5 

0.065 

0.086 

0.107 

35 


64.  Equation  36e  is  said  by  Takemitsu  to  improve  the  k  -  e  formulation 
performance  near  walls  since  it  includes  the  effect  of  diffusion  on  «  .  As 
evidence  of  its  superiority,  he  shows  that  in  a  trial  application  to  open 
channel  flow  (using  the  probable  constant  values  from  the  tabulation) ,  the 
resulting  eddy  viscosity  exhibits  a  minimum  at  the  bed  and  surface,  plus  a 
^hird  minimum  at  middepth,  a  reasonable  distribution  that  is  not  obtained 
with  some  approaches.  Equation  36e  does  not  include  a  bouyancy  term. 

Level  2-1/2  model 

65.  An  alternate  approach  to  turbulence  transport  modeling  is  presented 
by  Mel  lor  and  Yamada  (1982)  and  endorsed  for  numerical  flow  models  by  Blumberg 
and  Kantha  (1986).  The  level  2-1/2  model  describes  turbulent  kinetic  energy 
by 
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where  Sq  is  a  constant  defined  in  paragraph  66. 
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where 
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and  SH  is  defined  in  Equation  37i,  and  e  is  again  dissipation  but  is  a 
rate  expressed  as 


36 


« 


(37e) 


(2k)3/2 

*1 


where  Ax  is  a  length  scale. 

66.  The  level  2-1/2  model  also  includes  an  equation  for  the  master 
mixing  length  scale  i  ,  to  which  all  other  length  scales  are  proportional. 
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where 

p  -  density 

k  -  von  Karman  constant 
b  -  channel  width 

and  the  constants  and  their  values  are  listed  as  follows: 

Ax  -  0.92 

A2  -  0.74 

Bt  -  16.6 
B2  -  10.0 

Cj  -  0.08 

Ej  -  1.8 

E2  -  1.33 

Sq  -  0.2 

s,  -  0.2 

Level  2  model 

67.  The  complexity  of  the  level  2-1/2  model  is  daunting;  however,  the 
complexity  is  a  problem  only  to  the  extent  that  the  computer  computations  may 
be  too  slow  (or  too  expensive)  for  the  purpose.  A  significantly  simplified 
level  2  model  described  by  Mellor  and  Yamada  (1982)  may  be  appropri  e  for 
many  applications.  The  level  2  model  begins  with  the  level  2-1/2  model  and 
then  neglects  the  derivative  and  diffusion  of  turbulent  energy,  so  that 
Equation  37a  becomes 
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Equation  37f  for  the  length  scale  is  replaced  with  an  algebraic  expression 


(371) 


KXj  + 


where 


-  a 


m 

|  IXj!  V^k”  d  Xj 

«  — 
Jv^iTd  Xj 


(37m) 


38 


where  a  is  a  constant  approximately  equal  to  0.01  and  Equations  37h  and  3 7 i 
are  replaced  with 
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and  Rf  is  a  flux  Richardson  number 
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68.  The  attractiveness  of  the  k  -  c  and  Meilor- Yamada  models  is  that 
they  incorporate  implicitly  many  of  the  processes  we  wish  to  capture  in  the 
calculation  of  boundary  layer  flows.  One  contributing  factor  that  they  do  not 
represent  is  the  effect  that  suspended  sediment  may  have  on  turbulence  damping 
apart  from  the  density  effect. 


Suspended  Sediment  Stratification 

69.  Coleman  (1981)  derived  an  expression  for  u(y)  in  uniform, 
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suspended  sediment,  stratified,  open  channel  flow  using  the  vertical  velocity 
profile  given  by  Scottron  (1967)  (notation  returns  to  the  x,y,z  system): 


u 

u. 


(38) 


where 


A  -  integration  constant 

Au/u.  -  channel  roughness  velocity  reduction  function  (Hinze  1959) 

(II//c)w(y/6)  -  wake  region  velocity  augmentation  function  (Coles  1956), 
where  w  is  a  function  of  y/6  ,  and  II  is  the  wake 
strength  coefficient 

In  Equation  38  v'  is  the  kinematic  viscosity  of  a  solid- liquid  two-phase 
mixture  at  an  elevation  y  above  the  bed,  given  by 


,  _  (1  -t-  2.5c  -*•  6.25c2  ♦  15.62c3) 
P  +  (pt  -  p)c 


(39) 


where  c  is  the  sediment  volume  concentration  and  p%  is  the  sediment  den¬ 
sity.  Thus,  the  effect  of  suspended  sediment  on  the  boundary  layer  structure 
is  explicit  in  Equation  38. 

70.  The  first  expression  in  parentheses  in  Equation  38  is  the  von 
Karman-Prandtl  law  of  the  wall,  while  the  last  two  terms  have  been  added  as  a 
result  of  subsequent  developments  In  boundary  layer  theory  and  constitute  the 
law  of  the  wake.  Coles  (1956)  stated  that  this  latter  law  is  characterized  by 
a  velocity  profile  at  a  point  of  separation  (or  reattachment),  at  which  the 
flow  is  locally  a  pure  wake  flow.  The  wake  region  in  turbulent,  fully 
developed  boundary  layers  is  an  indication  of  a  large-scale  eddy  structure 
that  is  constrained  primarily  by  inertia  rather  than  by  viscosity  (as  is  flow 
in  the  viscous  sublayer),  ~nd  therefore  is  similar  to  the  flow  region  in  a 
wake.  Coles  (1956)  found  the  following  empirical  relationship  for  w  : 


w 


2  sin2 


(90) 
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Evaluating  Equation  33  at  y/h  -  1  gives  the  expression  for  u^  ,  from 
which  the  velocity  defect  law  is  found  to  be 


>w-u 


u. 


(41) 


Coleman  stated  that  this  equation  describes  the  velocity  distribution  for  the 
entire  flow  depth  except  for  the  viscous  sublayer.  The  terms  inside  the 
square  brackets  are  the  logarithmic  part  of  the  von  Karman- Prandtl  velocity 
defect  law  and  an  intercept  term,  respectively.  The  value  of  the  wake 
strength  coefficient  II  can  be  determined  by  applying  a  straight-line  fit  to 
the  measured  data  points  from  the  asymptotically  logarithmic  part  of  the 
profile  to  the  high  values  of  y/h  ,  giving 


n 


K 

1 


U 


Jy/h-i 


(42) 


71.  Coleman  (1981)  found  an  empirical  relationship  between  II  and  the 
following  form  of  a  gross  flow  Richardson  number  (Monin  and  Yaglora  1971): 

Ri  -  i-Ph~Pt>)  (43) 

Pmu. 

where  ph  ,  pb  ,  and  pm  are  the  densities  of  the  sediment-water  mixture  at 
the  water  surface,  channel  bed,  and  depth -averaged  value,  respectively.  Thus, 
for  the  velocity  defect  profile  given  by  Equation  36,  the  suspended  sediment 
distribution  (and  therefore  the  effect  of  sediment- induced  stratification)  is 
only  implicitly  accounted  for  through  the  relationship  between  II  and  R4  . 
Unfortunately,  Coleman  was  not  able  to  verify  this  empirical  relationship  that 
he  found  from  laboratory  tests  using  any  other  data  sets.  In  addition,  the 
validity  of  applying  the  empirical  relationship  found  for  u  (given  by 
Equation  40)  in  all  bounded  shear  flows  has  not  been  yet  ascertained. 

72.  With  these  limitations  in  mind,  to  use  Coleman's  method  for 
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determining  u(y)  and  u.  ,  given  a  known  suspended  sediment  profile  c(y)  , 
Equations  40  and  41  would  have  to  be  iteratively  solved  along  with  a  digitized 
form  of  the  relationship  between  Rt  and  II  . 

73.  Ordonez  (1970)  applied  an  empirical  correction  to  the  velocity 
defect  law  to  account  for  deviations  (particularly  in  the  near-bed  zone)  of 
measured  velocity  profiles  in  sediment  stratified  flows  from  that  given  by  the 
von  Karman-Prandtl  velocity  defect  law.  His  modified  velocity  distribution  is 


Umax  -  u 


U. 


(I'*1"  i 


(44) 


where  is  an  empirically  determined  coefficient,  the  value  of  which  was 
found  to  vary  between  0.001  and  0.016.  Ordonez  found  a  correlation  between  ♦ 
and  k  ,  with  the  latter  given  by 


K 


1  +  ^(S  -  1) 

■TTT3c0 


(45) 


where 

cm  -  depth- averaged  sediment  concentration  by  volume 
cQ  -  maximum  (near -bed)  sediment  volume  concentration 
S  -  sediment  specific  gravity 

The  relationship  between  k  and  if  is  a  function  of  Ws/u„  ,  where  Ws  is 
the  sediment  fall  velocity.  Thus  to  determine  u(y)  and  u*  ,  Equations  44 
and  45,  along  with  the  implicit  relationship  between  k  ,  if  ,  and  u,  ,  have 
to  be  solved  iteratively. 

74.  Gust  (1976)  and  Gust  and  Walger  (1976)  reported  that  measurements  of 
mean  velocity  and  suspended  cohesive  sediment  profiles  in  a  tidal  channel  in 
the  North  Sea  yielded  deviations  from  the  law  of  the  wall  in  a  dilute 
seawater/clay  mineral  suspension.  The  same  result  was  obtained  from  labora¬ 
tory  experiments  in  a  recirculating  flume  in  which  tidal  flow  of  a  seawater/ 
clay  mixture  was  simulated.  Specifically,  deviations  from  the  logarithmic  law 
of  the  wall  were  found  for  10  <  y'  <  30,  where  y*  -  yu^/v  .  Gust  concluded 
that  the  law  of  the  wall  cannot  be  used  to  compute  bed  shear  stress  in  clay 
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suspension  flows.  Yet  he  did  not  present  an  alternate  method  for  determining 
the  velocity  profile,  and  therefore  u„  ,  in  flows  with  clay  suspensions. 

75.  Gust  also  found  that  the  measured  velocity  profiles  were  very  simi¬ 
lar  to  those  for  drag-reducing  polymer  flows.  Thus,  he  further  concluded  that 
drag  reduction  was  occurring  in  the  clay  suspension.  He  hypothesized  that 
drag  reduction  was  at  least  partially  caused  by  aggregates  of  clay  particles, 
which  change  shape  under  the  action  of  shear  or  vorticity  in  a  turbulent  flow. 

76.  Gust  and  Walger  noted  the  following  differences  between  clear  water 

(Newtonian)  flows  and  clay  suspension  (non-Newtonian)  flows:  (a)  the  viscous 
sublayer  thickness  for  clay  suspensions  was  about  five  times  larger  than  that 
for  clear  water  flows  under  the  same  flow  conditions;  (b)  the  values  of  k 
obtained  from  the  slopes  of  measured  clay  suspension  velocity  profiles  varied 
between  0.3  to  0.4;  (c)  the  velocities  in  the  logarithmic  layer  were  the  same 
for  the  two  types  of  flow  for  "comparable  Reynolds  numbers";  and  (d)  u,_ 

values  for  clay  suspension  flows  were  approximately  40  percent  less  than  those 
obtained  for  clear  water  flows . 

77.  With  regard  to  (b) ,  Gust  considered  the  variation  in  k  a  result  of 
the  mean  flow  experimental  error  of  5  percent,  and  not  due  to  the  presence  of 
suspended  cohesive  sediment.  Thus,  he  assumed  k  -  0.40  in  his  analysis. 
However,  many  researchers  have  found  that  k  is  a  function  of  the  rate  of 
turbulent  energy  dissipation,  which  in  turn  is  partially  governed  by  the  sus¬ 
pended  sediment  concentration  (Einstein  and  Chien  1954,  1955;  Einstein  and 
Abdel-Aal  1972;  Hino  1963;  Ippen  1971;  Wang  1981).  Coleman  (1981)  supported 
the  constancy  in  the  value  of  k  ,  and  attributed  the  conclusions  arrived  at 
in  these  studies  referenced  to  the  incorrect  manner  in  which  the  data  of 
Vanoni  (1946),  Einstein  and  Chien  (1954,  1955)  and  others  were  analyzed. 
Specifically,  Coleman  stated  that  the  existence  of  the  wake  flow  region  and 
its  effect  on  the  velocity  profile  were  unknown  at  that  time,  which  led  them 
to  erroneously  believe  that  it  did  not  matter  to  which  section  of  a  u(y) 
versus  log  (y)  plot  a  straight  line  was  fit  to  determine  the  value  of  k  . 

78.  Fortier  and  Scobey  (1926)  found  that  the  channel  velocities  for 
flows  carrying  suspended  cohesive  sediment  were  greater  than  for  clear  water 
flows.  Gust  (1976)  stated  that  this  observation  can  be  explained  in  terms  of 
the  occurrence  of  drag  reduction  as  follows:  the  thicker  viscous  sublayer, 
which  occurs  in  a  clay  suspension,  means  that  the  velocities  in  the  logarith¬ 
mic  layer  above  the  viscous  sublayer  must  be  greater  than  those  in  the  same 
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region  in  clear  water  flows  to  obtain  the  same  critical  u.  values  necessary 
for  erosion  of  clay  particles  from  the  bed  surface,  as  found  in  clear  water 
flows . 

79.  Gust  and  Southard  (1983)  studied  the  structure  of  two-dimensional 
equilibrium  boundary  layer  flows  using  two  types  of  bottoms  in  a  laboratory 
flume:  a  movable  smooth  bed  of  0.16 -mm  quartz  sand  and  a  nonmovable  concrete 
bottom  with  the  same  equivalent  sand  roughness.  They  found  that  the  measured 
velocity  profiles  did  not  follow  the  universal  law  of  the  wall  for  the  flows 
over  the  movable  bed,  in  which  sand  grains  were  rolling  over  a  smooth  sand  bed 
(weak  bed  load),  while  they  did  for  the  nonmovable  bed.  Further  comparison  of 
the  two  flows  revealed  that  for  the  flows  over  the  movable  sand  bed,  (a)  the 
logarithmic  layer  extended  further  into  the  wake  region;  (b)  the  flows  had  a 
reduced  value  of  k  -  0.32  +  0.04  (even  though  the  water  remained  sediment 
free);  (c)  the  logarithmic  layer  extended  down  to  the  top  of  the  moving  sand 
particles;  and  (d)  the  friction  diagram  showed  drag  reduction  that  was 
Reynolds -number  dependent. 

80.  With  regard  to  the  third  point,  Gust  and  Southard  found  that  with 
increasing  Reynolds  number,  the  viscous  sublayer  becomes  increasingly  weak¬ 
ened,  and  in  fact  is  eventually  replaced  by  a  layer  that  may  be  either  a 
"buffer  layer,  a  wake,  or  even  an  extension  of  the  overlying  logarithmic 
layer."  With  regard  to  the  first  and  third  points  collectively,  they  mention 
that  as  the  number  of  sand  grains  in  motion  increase,  the  logarithmic  layer 
extends  both  of  its  boundaries,  with  the  top  one  approaching  the  water  sur¬ 
face.  With  regard  to  the  fourth  point.  Gust  and  Southard  conjecture  that  drag 
reduction  is  controlled  mainly  by  the  bed  load.  In  conclusion,  Gust  and 
Southard  hypothesize  that  a  new  type  of  wall-bounded  shear  flow  theory  is  re¬ 
quired,  one  with  possibly  different  momentum  transfer  processes  and  a  velocity 
defect  region  that  extends  from  the  top  of  the  moving  sand  particles  to  the 
free  surface. 


Solution  Methods 

81.  The  boundary  layer  calculations  discussed  in  this  Part  range  from 
simple  but  crude  to  sophisticated  but  difficult.  Equations  9-13  are  based  on 
classical  steady- flow  analyses,  and  they  can  be  solved  analytically  if  the 
free -stream  velocity  and  boundary  roughness  are  known.  This  approach  has  been 
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commonly  used  in  field  and  laboratory  studies  of  estuarine  flows,  but  is  known 
to  be  deficient  under  many  circumstances  and  seriously  in  error  in  some  (e.g. , 
stratification  or  near  times  of  slack  water) . 

82.  A  more  rigorous,  but  still  analytical  solution  can  be  obtained  by 
using  a  multilayer  model  like  Kajiura's  (Equations  18-28)  for  unsteady  flows 
without  complications  of  density  stratification  or  sediment  transport.  For 
real-world  solutions,  with  these  complicating  factors,  it  is  not  presently 
clear  whether  the  additional  rigor  of  Kajiura's  solution  is  to  be  preferred 
over  specialized  modifications  to  the  steady- flow  equations  which  account  for 
stratification  (Equations  32-35). 

83.  Where  computational  resources  are  sufficient  and  the  necessary 
boundary  condition  data  are  available,  a  higher  order  turbulence  closure  is 
recommended.  Both  the  k  -  e  model  and  Mellor-Yamada  model  are  physically 
rigorous,  and  have  much  to  recommend  them,  but  the  authors  favor  the  levels  2 
and  2-1/2  Mellor-Yamada  for  estuarine  boundary  flows.  Their  use  requires 
numerical  solution  of  the  equations  using  approximate  methods  such  as  finite 
differences  or  finite  elements.  This  necessitates  use  of  a  computational  mesh 
with  appropriately  posed  boundary  conditions.  From  numerical  experiments  re¬ 
ported  in  the  literature,  it  seems  probable  that  extreme  care  is  required  in 
formulation  of  the  numerical  solution  (e.g.,  centering  of  differences). 
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PART  IV:  BOUNDARY  LAYERS  AND  STRESSES  DUE  TO  WAVES  AND  CURRENTS 


84.  Most  practical  estuarine  applications  involve  determination  of 
shear  stresses  due  to  both  waves  and  currents.  In  most  estuaries  wind¬ 
generated  surface  waves  would  fall  into  either  the  shallow-water  or 
intermediate -water  wave  classification.  Thus,  as  a  wave  propagates  across  the 
water  surface,  the  wave -induced  orbital  velocity  of  water  particles  near  the 
bed  induces  a  shear  stress  on  the  bed  surface  in  addition  to  the  shear  caused 
by  quasi-steady  tidal  currents.  Reasonable,  first-order  approximations  to  the 
orbital  velocities  and  particle  excursions  may  be  obtained  from  small  ampli¬ 
tude  (linear)  wave  theory.  Effects/processes  not  represented  by  linear  wave 
theory  include  wave  form  asymmetry,  net  mass  transport  in  the  direction  of 
propagation,  flow  through  porous  beds,  and  wave  breaking.  A  brief  review  of 
wave  mechanics  (mainly  small  amplitude  theory)  is  given  next. 

Water  Wave  Mechanics 


85.  For  a  small  amplitude  gravity  wave  of  period  T  ,  wave  length  L  , 
and  amplitude  a  ,  propagating  in  a  quiescent  water  body  of  mean  water  depth 
h  ,  the  wave  length,  wave  number  k  ,  and  circular  wave  frequency  a  are 
defined  by 


L 


tanh  (kh) 


(46) 


k-  2* 

T 


(47) 


a 


2w 

T 


(48) 


where  g  is  the  acceleration  of  gravity.  The  most  restrictive  assumption 
made  in  linear  wave  theory  is  that  the  wave  amplitude  is  "small"  relative  to 
the  wave  length  and  water  depth.  This  allows  linearization  of  the  kinematic 
and  dynamic  surface  boundary  conditions,  which  allows  the  boundary  conditions 
at  y  «  r)  +  h  ,  where  y  is  the  vertical  distance  coordinate  and  r?  is  the 
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displacement  of  the  water  surface  relative  to  the  mean  or  quiescent  water 
level,  to  be  applied  at  y  -  h  (quiescent  position  of  the  water  surface).  To 
quantify  "small,"  it  is  required  that 


ak  «  tanh  (kh) 


(49) 


86.  The  horizontal  velocity  component  of  water  particles  beneath  a 
single  monochromatic  wave  uQ  as  a  function  of  distance  above  the  bottom 
(i.e.,  y  -  0)  is  given  by 


agk  cosh  (ky) 
a  cosh  (kh) 


sin  (kx  -  at) 


(50) 


where  x  is  the  horizontal  distance  coordinate  and  t  is  time .  Thus ,  the 
velocity  is  predicted  to  vary  sinusoidally  in  time  using  small  amplitude 
theory.  This  is  not  a  valid  assumption  when  more  than  one  monochromatic  wave 
field  is  present  and  in  shallow  water  where  higher  harmonics  become  increas¬ 
ingly  dominant. 

87.  Water  particle  motion  beneath  a  symmetrical  (linear)  surface  wave 
follows  a  closed  elliptical  path,  with  the  amplitude  (one-half  the  particle 
excursion)  of  the  major  axis  (which  is  in  the  direction  of  wave  propagation, 
i.e.,  x-direction)  given  as 


Ao 


cosh  (ky) 
sinh  (kh) 


(51) 


Both  uQ  and  A,,  ,  the  horizontal  amplitude,  become  negligible  at  depths 
greater  than  about  one -half  a  wave  length. 

88.  Shear  stresses  exerted  on  the  seabed  by  short-period  (i.e.,  wind¬ 
generated)  waves  have  been  the  subject  of  substantial  research.  They  are 
usually  expressed  in  the  form  of  Equation  9,  with  Uj  ,  the  velocity  in  the  x- 
direction  at  a  specified  distance  above  the  bed,  being  the  maximum  wave  orbi¬ 
tal  velocity  near  the  bed  .  The  wave  friction  coefficient  f„  has  been 
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defined  by  Jonsson  (1965)  for  hydraulically  rough  turbulent  flow  to  be 


1 


♦  log 


1 

7 


-0.08  +  log 


Aob 


(52) 


where  is  the  orbit  amplitude  at  the  bed  (i.e..  Equation  51  evaluated  at 

y  -  h  );  and  k8  is  the  equivalent  bed  roughness  size.  Kamphuis  (1975) 
determined  the  following  explicit  approximation  for  fw  : 


0.4 


-0.75 


(53) 


which  agrees  reasonably  well  with  Equation  52  for  2  <  (A^/kg)  <  100  . 

89.  Linear  wave  theory  is  used  in  most  of  the  methods  presented  in  this 
Part  to  estimate  the  combined  bed  shear  resulting  from  the  passage  of  a  mono¬ 
chromatic  nonbreaking  wave  and  tidal  currents.  It  is  pertinent  to  point  out 
here  that  simply  adding  the  shear  stress  due  to  waves  to  that  due  to  currents 
is  not  adequate  for  most  applications,  for  the  same  reason  that  the  sum  of 
squares  is  not  equal  to  the  square  of  sums. 

Laminar  Flows 

90.  Collins  (1964)  derived  the  following  expression  for  the  near-bed 
velocity  profile  in  which  he  assumed  that  a  steady,  collinear  laminar  current 
was  not  affected  by  oscillatory  wave-induced  motion: 


u  *  U*  [cos(kx  -  at)  -e^  cos(kx  +  0y  -  0t)  ]  ±  U  [0  A0y  -  0.04(/fy)2}  (54) 

where 

u  -  velocity  in  the  direction  of  the  x-axis 

U«  -  amplitude  of  the  horizontal  velocity  component  just  above  the  wave 
boundary  layer 

0  -  {o/2v)0i  where  v  is  the  kinematic  viscosity  of  water 
U  -  laminar  current  velocity  at  y  -  5/(<7/2i/)°  5 
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As  pointed  out  by  Sleath  (1984),  this  equation  is  a  first-order  est.mate  at 
best  since  it  linearly  superposes  the  laminar  flow  and  wave  shear  stress 
fields . 


Turbulent  Flows 

91.  This  section  discusses  bed  shear  stresses  and  boundary  layer  struc¬ 
ture  in  combined  low  frequency  (tidal)  and  high  frequency  (wind-generated  sur¬ 
face  waves)  turbulent  wave -induced  flows.  Three  methods  for  determining  bed 
shear  stress  and  boundary  layer  flow  in  turbulent,  combined  w i/e  and  current 
flows  are  discussed:  wave-current  friction  factors,  eddy  viscosity 
approaches,  and  turbulence  modeling. 

Bed  shear- friction  factor  approaches 

92.  Bijker  (1966)  assumed  a  log  velocity  profile  for  a  steady  current 
at  an  angle  to  the  direction  of  wave  propagation  and  developed  an  expression 
for  the  resultant  shear  of 


where 

rwe  -  resultant  combined  shear  due  to  waves  and  currents 
rc  -  shear  due  to  current  only 
p'  **  constant  coefficient  (given  as  0.45) 
k  -  von  Karman's  turbulence  constant 
C  -  Chezy  coefficient 

U  -  mean  (depth- averaged)  current  speed 
Swart  (1974)  found  that  p'  was  rot  a  constant,  but  was  the  following 
function  of  fw  : 


P' 


(56) 


Replacing  p'  in  Equation  55  with  Equation  56,  and  the  Chezy  coefficient  with 
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the  equivalent  current  friction  factor,  fc  ,  from  Equation  10  and  rearranging 
yields : 


rwc  -  0.5fcpU2  +  0.25f>4 


(57) 


where  p  is  the  density  of  water.  Bijker  used  waves  with  incidence  angles  of 
75  and  90  deg  and  a  steady  current.  Since  it  has  been  found  that  the  log 
velocity  profile  may  not  be  a  good  description  of  tidal  currents,  the  use  of 
the  log  law  in  Bijker 's  analysis  as  weir  as  the  use  of  steady  currents  in  his 
experiments  potentially  (at  least)  limits  the  usefulness  of  his  results. 

93.  Jonsson  (1966)  examined  bed  shear  due  to  combined  waves  and  cur¬ 
rents  over  plane  beds.  He  used  a  small  collinear  current  (U  <  u0>  and  an 
energy  analysis  to  find  that  the  instantaneous  value  of  shear  stress  is  equal 
to 


Twc  “  0.5fHCp(U  +  Uq)  |u+uj  (58) 

where  fwc  ,  the  friction  factor  for  combined  waves  and  currents,  is  given 
by: 


f  +  fCU 


(59) 


with  f„  given  by  Equation  52. 

94.  Experiments  by  Brevik  and  Aas  (1980)  showed  good  agreement  with 
Equation  59  for  0  <  U/u^,,  <  0.4  .  For  negative  values  of  U/uOT  (i.e.,  waves 
propagating  in  the  opposite  direction  of  the  current),  fwc  was  nearly  con¬ 
stant  at  about  0.35.  For  values  of  U/u^  greater  than  0.4,  the  computed 
values  of  fwc  were  somewhat  lower  than  that  given  by  Equation  59. 

95.  Experiments  used  by  Jonsson  for  verifying  Equation  58  were  for 
waves  propagating  in  the  direction  of  the  current  and  for  U  <  uD  ,  as  stated 
in  paragraph  93.  However,  as  Equation  58  does  not  explicitly  require  the 
second  assumption,  it  can  be  extended  to  the  general  case  (oblique 
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wave-currenc  interaction)  by  replacing  the  squared  velocity  sum  with  the 
square  of  the  resultant  velocity  due  to  the  superposition  of  U  and  uQ 
Equation  50  evaluated  at  x  -  0  gives 


Uo  -  u^  sin  (at) 

(60) 

while  U  and 

u0  expressed  in  terms  of  x-  and  z-components  give 

U2  -  U2  >  U2 

(61) 

2  2  2 
+ 

(62) 

The  resultant 

velocity,  ur  ,  can  be  expressed  as 

Ur2  “  +  ^ 

(63a) 

-  (Uj,  ♦  U^x)2  +  (U«  +  ^zf 

(63b) 

-  U2  +  u^  .  2<U)tut)I  ♦  UlUo,) 

(63c) 

Equation  58  gives  an  instantaneous  value  of  r„c  ;  therefore,  substituting 
Equation  63  into  Equation  58  and  integrating  over  one  wave  period  yields  the 
average  wave-  and  current- induced  bed  shear  stress: 


T 


*»c 


1 

T 


|  rMC(t)  dt 


(64a) 


f 

WC 

TF 


U2t  +  ^  (0 . 5at 
a 


0. 25  sin  (2crt) 


(64b) 


-  0 . 5fwcp(U2  ♦  0.5v&) 


(64c) 
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Because  of  the  form  of  Equation  58,  Equation  64c  cannot  be  separated  into  two 
parts  as  for  Equation  63.  It  should  be  noted  that  there  is  some  evidence  that 
the  resultant  velocity  for  a  quasi-steady  current  in  the  presence  of  oscil¬ 
latory  currents  is  not  a  linear  vector  sum  as  given  in  Equation  63a  (George 
and  Sleath  1979).  However,  the  phenomenon  is  not  well  enough  understood  to 
attempt  an  alternate  approach  here. 

96.  Both  Bijker's  and  Jonsson's  methods  (Equations  57  and  59,  respec¬ 
tively)  for  calculating  the  surface  shear  stress  contain  the  friction  factor 
due  to  currents  only  fc  .  Unfortunately,  as  discussed  in  Part  III,  there  is 
not  a  single  favored  way  to  compute  a  time -dependent  form  of  fc  for  bed 
surface  shear.  If  near-bed  velocity  measurements  are  available,  and  one  is 
willing  to  assume  the  validity  of  the  log  profile,  then  time-dependent  values 
of  fc  can  be  calculated  from  the  appropriate  expression  in  Table  1.  With 
the  lack  of  any  better  evidence,  it  may  be  necessary  to  use  Equation  10  with  a 
grain  size  ks  to  obtain  fc  or  to  use  one  of  the  more  rigorous  approaches 
discussed  previously,  back-calculate  fc  ,  and  substitute. 

97.  O'Connor  and  Yoo  (1988)  modified  Bijker's  model  to  account  for  the 
presence  of  the  wave  boundary  layer  and  current/wave  interaction,  in  particu¬ 
lar  reduction  of  the  current  speed.  They  divided  the  wave  period  averaged 
resultant  bed  shear  into  a  component  parallel  to  the  current,  <rb>p  ,  and 
normal  to  the  current,  <rb>n  ,  with 


<Tb>P  -  yTb 


(65) 


and 


<rb>n  -  /3rb 


(66) 


where  <  >  is  the  wave  period  operator;  rb  is  the  current- induced  time  mean 
bed  shear  stress;  and  7  and  ft  are  given  by 

7  -  O^OSq0-75^1’25  sin  (20)  (67) 
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(68) 


f 


ft 


or 


1  +  [0.36  -  0.14  cos  (26) 


where  8  is  the  angle  which  the  orbital  wave  motion  makes  with  the  direction 
normal  to  the  current;  and  p  is  the  relative  velocity  of  maximum  wave  orbi¬ 
tal  velocity  to  the  current  at  y'  ,  given  by 


with  uto  as  the  maximum  wave  velocity  at  the  top  of  the  wave  boundary  layer; 
and  y'  as  the  hypothetical  elevation  above  the  bed  (proposed  by  Bijker) 
where  a  line  emanating  from  the  seabed  (y  -  h)  is  tangential  to  the  current - 
only  velocity  profile.  If  the  wave-  and  current -induced  velocity  profile  can 
be  approximated  as  varying  logarithmically,  then  the  level  y'  is  approxi¬ 
mately  equal  to  eyj  ,  where  e  is  the  base  of  natural  logarithms,  and  y1 
is  the  height  above  the  (theoretical)  bed  at  which  the  velocity  is  equal  to 
zero  (Bijker  1966) . 

98.  The  term  a  in  Equations  67  and  68  is  defined  as  the  current 
velocity  reduction  factor,  which  accounts  for  reduction  of  current  in  the  wave 
boundary  layer,  and  is  evaluated  as 


a 


/?  +  S 


(70) 


where  6  is  the  boundary  layer  thickness  with 


and 


5  -  +  [o.43 


0 . 2  cos 


q  -  0 . 1  cos  (26)  -1.9 


(71) 


(72) 
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Obviously,  Equation  70  has  to  be  solved  iteratively  with  Equations  68  and  71 
for  a  .  O'Connor  and  Yoo  (1988)  suggest  using  a  Newton -Raphson  method. 
Comparison  of  this  method  with  several  data  sets  yielded  reasonable  agreement. 

99.  The  wave  friction  factor  f„  and  current-only  friction  factor  fc 
used  in  Equations  69  and  70  were  evaluated  by  O'Connor  and  Yoo  using  the 
following  equations: 


fc  -  0.016  ^J"1/3  (73) 


-0.194 

4  -  exp 

5.213 

TT 

V  J 

-  6.67 

for  >  2.0 

(74) 

f*  -  0.12  for  4^  <  2.0  (75) 


Eddy  viscosity  approaches 

100.  Methods  for  determining  the  instantaneous  and  wave  period  averaged 
bed  shear  stress  or  friction  coefficient  for  aombined  wave  and  current  fields 
as  a  function  of  both  time- invariant  and  time -dependent  eddy  viscosities  are 
also  available.  Whenever  possible,  equations  that  approximate  the  actual 
expressions  for  bed  shear  or  wave-current  friction  coefficient  are  presented 
for  ease  of  use. 

101.  Grant  and  Madsen  (1979)  presented  a  method  for  determining  (a)  the 
bed  shear  stress  on  a  rough  bottom  under  the  combined  action  of  waves  and 
currents  using  a  wave-current  friction  coefficient,  and  (b)  the  wave-  and 
current- induced  velocity  profile  both  inside  and  outside  the  wave  boundary 
layer  using  a  time -invariant  two-layer  eddy  viscosity  model.  They  used  the 
two- layer  eddy  viscosity  model: 


’*y  KJ  for  y  >  <5W 
«y  Kwcl  for  y  <  5W 


(76) 
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where 


i/t  -  eddy  viscosity 
y  -  height  above  the  bed 

|u.c|  -  magnitude  of  the  time-averaged  current  shear  velocity 
-  wave  boundary  layer  thickness 

|u„wej  -  magnitude  of  maximum  shear  velocity  due  to  combined  wave  and 
current  motion 

Thus,  u*wc  incorporates  wave-current  interaction  effects.  The  magnitude  of 
the  maximum  wave-  and  current- induced  bed  shear  stress  is  the  usual 

function  of  ju**,.)  : 


lTb«l  -  P  lu^cl2  “  P<*  l^xnl2 


(77) 


where  u^,  is  given  by  the  modulus  of  Equation  57  with  y  -  h  ;  and  a  is 
given  by 


a 


|u.l: 

TOT 


lu.l 

ThTT 


cos  4> 


(78) 


where  |ua|  is  the  magnitude  of  steady  current  velocity  at  a  height  y  -  a 
above  the  bottom;  and  <j>  is  the  angle  made  by  U,  with  the  wave  propagation 
direction. 

102.  The  implicit  equation  derived  by  Grant  and  Madsen  for  the  wave- 
current  friction  factor  fwc  in  Equation  77  is: 


'  ’ 

1/2 

2 

1/2 

0.097 

K 

+  2 

0.097 

kb 

K 

T^bTl 

J 

~^r 

l^obl 

J 

jt3/ * 

vf 


(79) 
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where 


kb  -  characteristic  dimension  of  the  physical  bottom  roughness 

jAobj  -  I'10)/0  (excursion  amplitude) 

K  -  turbulent  kinetic  energy 

<PC  -  angle  of  the  average  shear  stress 

and  V2  is  given  by  the  following  for  small  currents  relative  to  wave  action: 


V,  -  (2*> 


in’  i.r 


(80) 


When  the  currents  are  not  small  compared  to  the  wave  motion,  the  general  form 
of  V2  (given  by  Equation  14  in  Grant  and  Madsen  (1979))  must  be  used.  The 
K  term  in  Equation  79  is  equal  to 


K  - 


1  _ 1 _ 

2FF2  (Ker2  2f  "2  ♦  Kei2  2f"2)1/2 


(81) 


where  -  kb/30i'  in  which  X'  -  (K  |u*  I  )/a  and  Ker  and  Kei  are 

WC 

Kelvin  functions  of  zeroth  order.  The  term  X'  is  a  measure  of  the  wave 
boundary  layer  thickness.  Once  Equation  79  has  been  solved  for  fwc  ,  Equa¬ 
tion  77  can  be  used  to  evaluate 

103.  Tanaka  and  Shuto  (1981)  used  a  one-layer  time- invariant  eddy  vis¬ 
cosity  model  to  derive  a  linear  equation  for  the  wave -current  friction 
coefficient.  The  eddy  viscosity  they  used  is  given  by 


v 


t 


"^•wc 


(82) 


104.  The  bottom  shear  stress  in  terms  of  the  eddy  viscosity  is  given  by 
the  familiar  expression: 


Pv  t 


3u 

T9 


pku^y 


3u 

Ty 


at  y  -  yD 


(83) 
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where  yc  is  the  roughness  height  and  u  -  uc  +  u„  ,  the  linear  combination 
of  the  boundary  layer  velocity  profile  due  to  currents  and  waves.  Equations 
for  uc  and  u*  are  given  by  Tanaka  and  Shuto.  They  derive  an  approximate 
(simplified)  implicit  equation  for  fwc  on  a  rough  bed  that  is  valid  for 
Ub./^yo  >  50  : 


=  2B  (i)2  *  4BC  (i)  C°S  *  +  °2 


(84) 


with 


B  = 


In 


-  1 


r  1-1/2 


c  =  *  4 

0.25  ♦  0.101 

In 

yc 

-0.5  In  (fwc)+2.42 

2 

n 

k 

Ub- 

J 

4 

(85) 


(86) 


105.  Tanaka  and  Shuto  also  derive  an  expression  for  fwc  that  is  valid 
for  wave-  and  current- induced  turbulent  flows  over  smooth  beds : 


f~  • 2p2  (y  * 4rc  [t]  c°s  *  *  gz 


(87) 


with 


F 


0.4 


--1 

/XT  *  0.852 


(88) 


' 

0.25  +  0.101 

2  In 

1 

-  In  +  n  Sfifl 

l 

vXc 

> 

[  *  ) 

4 

(89) 
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where  ^  -  2A<,  .  For  the  smooth  bed  formulation,  the  thickness  of  the  vis¬ 
cous  sublayer  was  used  in  place  of  the  roughness  height.  Measured  wave-  and 
current-induced  shear  stresses  in  a  laboratory  flume  showed  fairly  good  agree¬ 
ment  with  those  predicted  using  Equation  83. 

106.  Tanaka,  Chian,  and  Shuto  (1983)  extended  the  work  of  Tanaka  and 
Shuto  (1981)  by  using  a  time-invariant  two-layer  eddy  viscosity  model.  They 
used  the  following  two-layer  formulation  suggested  by  Mellor  and  Gibson 
(1966): 


xyu.wc  for  0  <  y  <  d 

/cdu.  for  d  <  y  <  h 

*«c  J 


(90) 


where  d  is  the  height  of  the  inner  layer.  Wave-current  interactions  were 
limited  to  opposing  currents  (^  of  180  deg)  and  favorable  currents  (0  of 
0  deg) . 

107.  There  are  two  basic  differences  between  this  work  and  that  of 

Grant  and  Madsen  (1979):  (a)  the  eddy  viscosity  as  given  by  Equation  90  is 

continuous  across  the  interface  between  the  inner  and  outer  layers;  and 

(b)  the  effect  of  a  steady  current  is  included  in  determining  the  height  of 
the  inner  layer,  whereas  it  was  not  explicitly  accounted  for  in  the  work  by 
Grant  and  Madsen, 

108.  The  expression  derived  by  Tanaka,  Chian,  and  Shuto  for  fwc  using 
the  two-layer  eddy  viscosity  model  is  very  complicated  compared  with  that  ob¬ 
tained  with  the  one-layer  model  (Equation  84),  and  yet  the  difference  between 
the  two  models  (for  resulting  fKC)  is  insignificant.  In  addition,  comparison 
of  predicted  velocity  profiles  using  the  one-layer  and  two-layer  models  with 
profiles  measured  in  a  wind  tunnel  showed  that  in  most  cases  profiles  pre¬ 
dicted  by  the  one-layer  model  were  closer  to  the  measured  profile  than  those 
obtained  with  the  two-layer  model. 

109.  Madsen  and  Wikramanayake  (1988)  presented  an  analytical  model  for 
turbulent  wave-current  bottom  boundary  layer  flow  using  a  three-layer 
time-invariant  eddy  viscosity  formulation.  They  started  by  linearizing  the 
two-dimensional  (vertical  plane)  boundary  layer  equations  for  an  incompres¬ 
sible  homogeneous  fluid  over  a  plane  boundary  using  an  order-of-magnitude 
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analysis.  Next,  by  expressing  u  (horizontal  flow  velocity)  and  p  (pres¬ 
sure)  as  the  sum  of  time-invariant  (current)  and  time-varying  (wave)  compo¬ 
nents,  they  separated  the  linearized  boundary  layer  equation  into  the  follow¬ 
ing  two  equations  for  small-amplitude,  monochromatic  progressive  waves  and  a 
steady  current,  respectively: 


It 


0  * 


(91) 


(92) 


where  vp  is  the  mean  pressure  gradient. 

110.  Madsen  and  Wikramanayake  then  assumed  a  complex  harmonic  function 
solution  for  the  linear  wave  boundary  layer  equation  (Equation  91)  and  derived 
the  following  equation  for  the  wave-induced  orbital  velocity  uj  within  the 
wave  boundary  layer: 


d 

Hy 


(‘'t  -  i^u,,  -  0 


(93) 


where  u^  is  a  complex  function  of  y  and  i  is  yJ-X 

111.  Integration  of  Equation  93  from  the  bottom  (y  “  0) ,  where 
i>t(3u/dy)  -  rc/p  ,  to  a  level  y  yields 


•'t 


(94) 


The  pressure  gradient  term  on  the  right-hand  side  of  this  equation  is  then 
dropped  as  it  is  small  relative  to  rc/p  ,  which  results  in  the  law  of  the 
wall  for  the  equation  governing  the  current  boundary  layer: 
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c 


(95) 


du„ 


"3y 


u. 


with  u*c  as  the  shear  velocity  associated  with  the  current  boundary  shear 
stress. 

112.  Madsen  and  Wikramanayake  (1988)  represented  the  time-invariant 
eddy  viscosity  term  included  in  Equations  93  and  95  using  the  following 
formulation: 


't  -  i 


*u.  y 

•wc-' 


L  *u.cy 


y  ^  atS 

atS 


(96) 


*u.  a, S  a,6  S  y  < 


a,  S 


*  y 


where 

a  -  fraction  of  the  wave  boundary  layer  thickness  over  which  the  eddy 
viscosity  varies  linearly 

S  -  scale  of  the  wave  boundary  layer  thickness 

e  -  u*c/u*wc  »  which  represents  the  relative  magnitude  of  the  current 

turbulence  intensity  to  the  combined  wave-current  turbulence 
intensity 

These  definitions  will  apply  to  the  discussion  of  Madsen  and  Wikramanayake ' s 
method  in  paragraphs  112-121.  The  advantages  of  this  three-layer  eddy  viscos¬ 
ity  model  over  the  two-layer  models  previously  discussed  are  that  (a)  it  is  a 
continuous  function  in  y  ,  and  (b)  it  is  applicable  for  both  weak  and  strong 
current-wave  interactions. 

113.  After  incorporating  Equation  96  into  Equation  93,  the  latter  is 
solved  for  the  three  domains  specified  in  Equation  96.  The  solution  for  the 
near-bottom  (y  <  a6)  domain  is 


ud 


+  i  Kei 


♦  B 


♦  i  bei 


(97) 
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where  Ker  ,  Kei  ,  ber  ,  and  bei  are  Kelvin  functions  of  zeroth  order; 
f  “  y/5  ,  with  6  -  Ku*  /<r  ;  and  A  and  B  are  arbitrary  complex  constants 
(to  be  specified  by  boundary  conditions  later). 

114.  For  the  intermediate  domain  (ai  <  y  <  a6/e),  the  solution  of 
Equation  93  is 

Uj  «  C  +  De'SUa!  (98) 


where  ^1  -  (1  +  i)/^2  ;  and  C  and  D  are  complex  constants. 

115.  For  the  outer  domain  (y  >  a S/e),  the  solution  of  Equation  93  is 

ud  -  4Ker(2^7r)  +  i  Keife^T*")]  +  fjberk/fTT)  ♦  i  beifeVfA")]  <9' 


where  E  and  F  are  complex  constants. 

116.  The  conditions  applied  to  Equations  97-99  are  the  no-slip 
conditions  at  y  -  y0  -  k./30  ,  the  approach  to  the  free-stream  velocity  as 
y  -*  <*>,  and  matching  conditions  between  the  three  domains.  Applying  these 
conditions  resulted  in  the  following  five  linear  equations: 

A  [Ker  (2^)  ♦  i  Kei  (2^)]  ♦  B  [berfc^)  ♦  i  bei  (2^)]  -  -u*  (100) 


A  Ker 


M  +  i  Kei  (2^)]  *  B  [ber  (2^0)  +  i  bei  <2^)1 


(101) 


^ia  ~^ia 

C  e  +  D  e 


i/la  ,l  ♦  D  e"V,i“  /e  =  E[Ker(2Va/€)  +  i  Kei(2yk/e)] 


(102) 
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AlKer'feyk)  +  Kei'  (2VS)]  •  4b«'(*£)  +  1  bei'  (2^)] 


_ _ _ *  V~ia 

C  eyfie*  -  D/Ii" 


(103) 


n — ^ ia  It  I ia  /t 

C  eyie  -  Dyie 


E  Ker 


(2^/e)  +  i  Kei'  [lifa/t] 


(104) 


where  ub  is  the  near-bottom  velocity.  These  five  equations  can  be  solved 
for  the  constants  A,  B,  C,  D,  and  E  once  the  following  five  parameters  are 
specified: 

K  ,  ub  ,  a  ,  u,c  ,  u.wc  (105) 


Madsen  and  Wikramanayake  assume  that  ka  ,  Uj,  ,  and  u*c  are  either  known  or 
specified,  which  leaves  u*wc  and  o  to  be  determined.  They  regard  a  as  a 
"free  (fitting)  parameter,"  and  uses  a  closure  hypothesis  to  determine  u, 
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117.  The  closure  hypothesis  Madsen  and  Wikramanayake  used  assumes  that 


u*wC  =  ^7p 


(106) 


with  t„  the  maximum  combined  wave-current  bed  shear  stress  vector  expressed 
by 


141  =( 


ri  t\  *  2 rcr, 


\l/2 

cos  *„c) 


*  r. 


(l  ♦  2a*2  cos  tfwc  +  /**) 


1/2 


(107) 


where 

r„  -  maximum  wave  bottom  shear  stress 

rc  -  average  current  bottom  shear  stress  at  an  angle  <£wc  to  the  wave 
direction 

H  -  u„e/u*w  ,  and  thus  equal  to  the  relative  magnitude  of  current  to 
wave  bottom  turbulence  intensity 
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Combining  Equation  107  and  the  definitions  for  e  and  n  yields 


«  =  (l  +  2/i2  cos  ^wc  ♦  p4)  1M 


(108) 


118.  Madsen  and  Wikramanayake  define  a  wave  friction  factor,  fwc  , 
using  Equation  10  with  r  -  rw  ,  and  -  Uj,  ,  and  a  modified  wave  friction 
factor,  f*c  ,  defined  to  be 


Co  -  fwc  f1  +  2^2  c°s  Cc  +  M4)  1/2 


(109) 


Madsen  and  Wikramanayake  obtain  the  following  implicit  relationship  for  the 
modified  friction  factor: 


|A[Ker'  (2^)  +  i  Kei'  (2/^)] 

“*>  (110) 

♦  B[ber'  (2/?7)  ♦  bei'  (2J?^)  ]  |  for  limit  f  *  f0 


in  which 


fcn/30 


V30 


+  2#i2  cos 
— 


♦  M*)* 


K,  i 


(in) 


where  k„  is  the  boundary  roughness  size 
and 


Ab  =  (l  +  2 A*2  cos  ^HC  +  ^4)1/Z  =  Ab  (l  +  2/i2  cos  ^wc  +  p4)  /Z  (112) 

119.  The  current  problem  (Equation  9?)  is  solved  and  the  no-slip  bottom 
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boundary  condition  and  matching  conditions  between  the  three  regions  are 
applied  to  obtain 


uc  =  —  —  In  2.  for  y  <  a6 

u.wc  *  y0 


(113) 


u 


c 


for  aS  <  y  <  aS/e 


for  y  >  a6/e 


(114) 


(115, 


Madsen  and  Wikramanayake  (1988)  point  out  that  the  solution  for  the  outer 
region  (Equation  115)  is  limited  by  the  validity  of  the  law  of  the  wall,  and 
thus  should  not  be  used  for  y  greater  than  a  fraction  of  the  current 
boundary  layer  thickness. 

120.  Madsen  and  Wikramanayake  (1988)  outline  an  iterative  technique  for 
solving  the  wave-current  problem.  The  procedure  involves  initially  assuming 
that  u.e  is  much  smaller  than  u*„  and  u^,c  .  Thus,  in  the  limit,  both  p 
and  e  are  equal  to  zero.  Equations  110-112  are  used  to  determine  the  value 
of  f'c  ,  and  then  fwe  from  Equation  109.  Using  this  value  of  fwc  ,  rw  is 
determined  using  Equation  10,  from  which  new  estimates  of  /j  and  e  are 
obtained  using  the  definition  for  p  and  Equation  108.  These  new  values  for 
H  and  e  then  replace  the  original  estimates  for  these  parameters  in  the 
iterative  procedure  described  previously.  Iteration  proceeds  until  conver¬ 
gence  is  achieved.  A  value  for  a  is  assumed  in  this  method.  Then  the  solu¬ 
tion  to  the  wave-current  problem  can  be  obtained  by  solving  Equations  100-104 
(with  the  calculated  values  for  n  ,  t  ,  u*wc  and  5)  for  the  wave  problem 
and  Equations  112-115  for  the  current  problem. 

121.  Madsen  and  Wikramanayake  (1988)  compared  the  results  from  their 
analytical  model  with  those  obtained  by  Davies,  Soulsby,  and  King  (1988)  (dis¬ 
cussed  in  the  next  section)  using  a  numerical  turbulence  closure  model. 
Excellent  comparison  was  achieved  with  a  value  of  a  of  1/2  .  This  is  par¬ 
ticularly  significant  considering  the  relative  simplicity  of  Madsen  and 
Wikramanayake' s  model  as  compared  with  the  full-scale  turbulence  closure  model 
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developed  by  Davies,  Soulsby,  and  King. 

122.  The  largest  difference  seen  between  the  results  of  Madsen  and 
Davies  was  that  the  former  showed  much  less  sensivity  to  the  angle  between  the 
direction  of  wave  propagation  and  the  current  than  the  latter.  For  this  rea¬ 
son,  Madsen  has  initiated  development  of  a  time-varying  eddy  viscosity  model 
(Madsen  and  Wikramanayake  1988).  This  model  has  not  yet  been  completed. 

123.  Lavelle  and  Mofjeld  (1983)  developed  a  flow-eddy  viscosity  inter¬ 
dependent  model  of  the  bottom  boundary  layer  using  a  time-varying  eddy 
viscosity.  While  their  model  does  not  include  effects  of  combined  wave  motion 
and  currents,  they  did  evaluate  the  model  for  the  case  of  oscillatory  flow  at 
tidal  frequencies.  They  found  that  if  time  variations  in  the  eddy  viscosity 
are  neglected,  then  the  maximum  bed  shear  stress  is  underestimated  and  the 
velocity  profile  is  distorted,  particularly  near  flow  reversal.  Some  aspects 
of  their  model  are  described  in  the  following  paragraphs. 

124.  The  time-variant  eddy  viscosity  used  had  the  form: 


*yK(t)  | 
kS |ub(t) | 


for  yc  <  y  <  5 
for  y  >  6 


where 

6  =“  log  layer  thickness 

uj,(t)  -  magnitude  of  modified  shear  velocity,  given  by: 


K(t)  |  =  [u?(t)  ♦  «/t  u*(t  *  T/4) j1/2 


(116) 


(117) 


where  u*  is  the  shear  velocity  defined  by: 


u.  = 


rb 


(118) 


The  bed  shear  rb  as  a  function  of  time  is  evaluated  using  the  first  expres¬ 
sion  on  the  right-hand  side  of  Equation  83. 
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125.  One  result  from  Lavelle  and  Mofjeld's  model  was  that  rb  and 
U0  ,  the  free-stream  velocity,  are  not  in  phase.  They  presented  the  following 
equation  for  rb  that  incorporates  the  relative  phasing  of  rb  and  UQ  ,  and 
both  linear  and  quadratic  friction  laws: 


rb(C/T)  -  g*!U0(t/T  *  6)  |*  Uc(t/T  ♦  9)  (119) 

where  9  is  the  fraction  of  the  wave  (tide)  cycle  between  the  zeros  of  rb 
and  U0  ;  P  has  a  value  between  0  and  1  (e.g.,  when  p  -  1 ,  rb  is 
quadratic  in  U„) ;  and  is  a  coefficient  (e  g. ,  gtf  is  analogous  to  CD 

(drag  coefficient)  when  p  -  1).  The  bottom  friction  was  also  found  to  be 
(implicitly)  dependent  on  the  degree  of  temporal  variation  of  the  eddy 
viscosity. 

126.  Lavelle  and  Mofjeld  also  found  that  flow  acceleration  and 
deceleration  cause  the  steady  state  near-bed  log  velocity  profile  to  be  modi¬ 
fied.  Their  theoretical  near-bed  velocity  function  contains  additional  terms 
in  y  In  (y/yQ)  and  y  .  Soulsby  and  Dyer  (1981)  confirm  this  finding.  This 
emphasizes  the  previously  stated  finding  that  use  of  conventional  log  velocity 
profiles  to  fit  accelerating  flow  measurements  may  result  in  inaccurate  esti¬ 
mates  of  the  shear  velocity  and  bottom  roughness.  This  also  points  out  possi¬ 
ble  limitations  of  the  previous  methods  reviewed  in  this  section,  since  most 
use  a  log  velocity  profile  for  the  steady  current. 

Turbulence  modeling  approaches 

127.  The  authors  recommend  the  five  papers  on  "Turbulence  Modeling  of 
Surface  Water  Flow  and  Transport"  which  appear  in  the  August  1988  American 
Society  of  Civil  Engineers  Journal  of  Hydraulic  Engineering.  These  papers 
review  the  state  of  turbulence  modeling,  and  thus  would  greatly  benefit  the 
reader  in  assessing  the  material  discussed  herein.  Another  recommended  paper 
is  that  of  Mellor  and  Yamada  (1982),  which  was  discussed  in  Part  III.  As  it 
is  beyond  the  scope  of  this  report  to  summarize  all  the  material  included  in 
even  these  six  papers  mentioned,  the  discussion  that  follows  is  limited  to  a 
description  of  one  turbulence  closure  model  that  has  been  used  to  specifically 
model  nonstratified  wave-  and  current-induced  boundary  layer  flow.  The  model 
reviewed  is  not  the  most  complete  or  advanced,  certainly  not  the  most  numeri¬ 
cally  sophisticated;  nevertheless,  it  is  able  to  reproduce  known  oscillatory 


boundary  layer  structure  using  a  seemingly  far  less  complicated  scheme.  As 
such,  it  appears  to  be  particularly  attractive  for  inclusion  in  a  full-scale 
three-dimensional  turbulence  model.  Two  closure  models  that  account  for  sedi¬ 
ment  stratification  are  reviewed  in  the  next  section. 

128.  Davies,  Soulsby,  and  King  (1988)  describe  a  relatively  simple 
turbulence  closure  model  of  oscillatory  rough-turbulent  boundary  layer  flow 
that  is  used  to  model  wave-current  interaction  in  the  near-bed  zone.  The 
philosophy  used  by  Davies,  Soulsby,  and  King  was  "to  model  the  flow  in  a  way 
no  more  complicated  than  can  be  warranted  by  present  abilities  to  make  con¬ 
trolled  field  measurements."  Advantages  and  features  of  their  approach  over 
similar  models  include  the  following:  (a)  it  can  be  applied  to  any  wave- 
current  condition,  from  large  waves  on  a  weak  current  to  small  waves  on  a 
strong  current;  and  (b)  the  turbulence  closure  scheme  involves  no  new  assump¬ 
tions  other  than  those  that  have  been  used  in  separately  modeling  wave  and 
current  boundary  layers. 

129.  Several  assumptions  are  incorporated  in  the  model:  (a)  the  flow 
is  horizontally  uniform  between  the  flat  rough  bed  and  horizontal  water  sur¬ 
face;  (b)  there  is  no  surface  shear  stress;  (c)  boundary  layer  approximation 
(k$wc  «  1)  is  valid;  and  (d)  waves  have  a  small  amplitude  and  satisfy  the 
Stokes  condition  (ka  «  1) .  This  last  assumption  allows  for  neglecting  the 
convective  acceleration  terms  in  the  momentum  equations. 

130.  The  following  procedure  is  used  in  the  model  to  initialize  bound¬ 
ary  layer  flow:  (a)  a  steady  uniform  starting  current  is  generated;  (b)  wave 
motion  of  a  given  frequency  and  incident  (with  respect  to  the  current  direc¬ 
tion)  angle  is  then  superposed  on  the  current;  and  (c)  the  model  is  run  to 
allow  convergence  to  the  combined  wave-current  steady  state.  The  justifica¬ 
tion  given  for  starting  with  a  steady  current  is  that  it  provides  consistent 
initial  vertical  profiles  for  velocity  and  turbulence  parameters. 

131.  The  linearized  horizontal  momentum  equation  in  the  z-direction 
(current  direction)  that  is  solved  to  generate  the  starting  current  is 

aw  _  _  i  ap  a  (120) 

It'  -pTz  Sy  [~J 


where  w(y,t)  is  the  mean  velocity  component  in  the  z-direction,  and  Tyz  is 
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the  Reynolds  shear  stress  in  the  z-direction  associated  with  the  current.  The 
latter  is  given  by 


y* 


Pvt 


dw 

'Sy 


(121) 


where  i/t(y,t)  is  determined,  as  described  in  the  following  paragraphs,  using 
a  turbulence  closure  approach. 

132.  The  following  two-dimensional  momentum  equations  are  used  to 
superpose  wave  motion  on  the  steady  current: 


du  _  l  dp  +  a  fO  (122) 

It  P  Tx  Fy  |7J 


aw  =  .  i  ap  a  (123) 

St  p  oz  ly  [  p  J 


where  r ^  and  ryz  are  components  of  the  shear  stress  in  the  x  and  z 
directions,  given  by 


y* 


P»  t 


3u 

Sy 


y* 


Pv  t 


3w 

Sy 


(124) 


133.  The  pressure  gradients  in  Equations  122  and  123  are  chosen  to 
yield  the  wave  motion  specified  in  the  following  equations  in  the  region  above 
the  wave-current  boundary  layer  (i.e.,  y  >  6WC)  : 


1  dp  m  dUx 

p  Sx  dt 


i  ap  =  duz 

p  Sz  dt 


+  Vp 


with  Ux  =  U0  sin  4>  cos  (at) 


with  Ut  =  U0  cos  4>  cos  (at) 


(125) 


(126) 
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where 


UItUt  -  horizontal  water  particle  velocity  components 

134.  In  the  solution  to  Equations  122-126  (together  with  appropriate 
boundary  conditions),  to  add  waves  to  the  already  established  steady  current, 
vp  was  set  equal  to  a  constant.  The  model  was  run  until  the  new  wave-current 
equilibrium  was  achieved.  Davies,  Soulsby,  and  King  (1988)  mention  that  set¬ 
ting  Vp  equal  to  a  constant  is  not  consistent  with  the  assumed  horizontal 
water  surface.  However,  they  add  that  the  pressure  gradients  of  concern  are 
usually  at  least  two  orders  of  magnitude  smaller  than  the  gradients  associated 
with  waves.  Therefore,  the  implied  surface  slope  is  negligible  and  can  be 
thought  of  as  an  artifact  of  the  numerical  scheme. 

135.  The  turbulence  closure  scheme  Davies,  Soulsby,  and  King  used  to 
solve  for  i/t(y,t)  is  described  in  detail  by  King,  Davies,  and  Soulsby 
(1985).  The  equation  for  turbulent  energy  is  given  by: 


k  _  M  fauf  fawf  .  .  a  iy  a k 

t  "Vt  pyj  M  Jy^Jy 


(127) 


where  K(y,t)  is  the  turbulent  energy  and  a'  is  the  ratio  of  eddy  diffusi- 
vities  of  turbulent  energy  K  and  momentum. 

136.  The  boundary  conditions  for  Equation  127  represent  zero  energy 
flux  across  the  water  surface  and  bed  surface.  The  turbulence  closure  scheme 
is  closed  using  turbulence  scaling  laws  for  vt  and  c  : 


i/t  =  c0fK1/2 


(128) 


where  cQ  and  Ci  are  constants  and  i  is  the  mixing  length,  given  by 


I  =  -*  K1'2  3  ft'2  A 

K  ~r  iy  [t  J 


(129) 


137.  For  the  steady  state  starting  current,  Equations  120,  121,  and 
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127-129  are  solved  together,  while  for  the  superposition  of  waves  and  cur¬ 
rents,  Equations  122-129  are  solved.  A  log  linear  depth  transform  was  applied 
to  these  equations,  after  which  Equations  120,  122,  123,  and  127  were  solved 
iteratively  using  an  implicit  Crank-Nicholson  scheme.  The  convergence 
criterion  was  based  on  the  predicted  vt  distribution.  Davies,  Soulsby,  and 
King  (1988)  used  the  values  given  by  Vager  and  Kagan  (1969)  for  the  constants 
in  Equations  120-129:  c0  -  c1/4  and  C!  -  c3/k  ,  with  c  -  0.046,  a  -  0.73, 
and  k  -  0.4.  King,  Davies,  and  Soulsby  (1985)  found  the  model  to  be  rather 
insensitive  to  changes  in  the  values  of  c  ,  a  ,  and  k. 

138.  The  model  was  used  to  simulate  oscillatory  boundary  layer  flow 
with  water  depth  of  10  m;  bed  roughness  of  0.5  cm;  steady  surface  current 
speed  of  1  m/sec;  8-sec  waves  with  incident  angles  of  0,  45,  and  90  deg;  and 
near-bottom  velocity  amplitudes  of  0.5,  1.0,  and  1.5  m/sec.  Distributions  of 
velocity,  turbulent  energy,  eddy  viscosity,  and  shear  stress  for  wave-current 
interaction  cases  were  presented.  Simulations  with  a  steady  current  and 
oblique  waves  showed  that  the  instantaneous  velocity  vector  varies  in  a  highly 
nonlinear  manner  with  the  amount  of  turbulent  energy  acting  to  retard  the 
near-bottom  flow.  It  was  also  found  that,  as  a  result  of  the  different 
amounts  of  turbulent  energy  generated  in  the  two  halves  of  the  wave  cycle,  the 
wave  cycle  averaged  horizontal  velocity  veers  away  from  the  direction  of  wave 
propagation.  Thus,  the  final  angle  between  the  current  and  waves  is  greater 
than  the  initial  oblique  angle.  The  degree  of  veering  varied  significantly 
over  the  water  depth.  For  example,  in  one  simulation  with  an  initial  incident 
wave  angle  of  45  deg,  the  veering  angle  of  the  mean  velocity  vector  varied 
from  15  deg  near  the  bottom  to  4  deg  at  the  water  surface .  The  depth-averaged 
current  veered  4.6  deg,  which  would  definitely  affect  sediment  transport  cal¬ 
culations  under  wave-current  combined  flows.  The  model  also  simulated  the 
nonlinear  amplification  of  the  bed  shear  stress  and  the  increase  in  the  oscil¬ 
latory  wave-current  boundary  layer  thickness  resulting  from  wave-current 
interaction.  Details  of  the  case  studies  are  given  by  Davies,  Soulsby,  and 
King  (1988). 

Stratification  effects 

139.  The  degree  to  which  salinity  and  suspended  sediment  stratification 
affect  the  wave-  and  current-induced  flow,  as  well  as  the  ability  of  the  flow 
to  stratify,  are  strongly  coupled  to  wave  motion  (Grant  and  Madsen  1986).  For 
example,  the  relatively  large  bed  shear  stresses  generated  by  wave-induced 


70 


oscillatory  motion  often  result  in  resuspension  of  large  quantities  of  sedi¬ 
ment,  which  may  cause  the  flow  to  stratify.  Nevertheless,  Glenn  and  Grant 
(1987)  found  that  sediment  stratification  does  not  occur  within  the  wave 
boundary  layer  for  typical  continental  shelf  conditions  due  to  small  concen¬ 
tration  gradients  usually  present.  However,  suspended  sediment  gradients 
under  estuarine  flow  conditions  may  be  one  to  two  orders  of  magnitude  greater 
Chan  those  that  occur  on  the  shelf,  and  thus  can  greatly  affect  the  wave-  and 
current- induced  boundary  layer  flow.  Efforts  at  modeling  the  effects  of 
salinity  and  suspended  sediment  stratification  on  bed  shear  and  boundary  layer 
flows  in  combined  waves  and  currents  are  reviewed  here. 

140.  Hagatun  and  Eidsvik  (1986)  describe  a  turbulence  closure  model 
that  simulates  sediment  stratified  oscillatory  turbulent  boundary  layer  flows. 
The  following  assumptions  are  incorporated  in  their  model:  negligible  inter¬ 
action  effects  between  suspended  sediment  particles,  except  for  the  fall 
velocity  component;  suspended  particles  advected  at  the  local  flow  velocity; 
collinear  waves  and  currents;  settling  velocity  of  sediment  particles  set 
equal  to  a  constant;  horizontal  scale  of  the  mean  flow  assumed  to  be  much 
greater  than  the  thickness  of  the  wave-current  boundary  layer,  so  that  hori¬ 
zontal  convective  acceleration  and  mean  vertical  velocity  terms  can  be 
neglected  in  the  momentum  equations.  The  equations  for  momentum  and  mass 
conservation  then  reduce  to 


3p u  _  _  3p  +  a 
at  3x  Jy 


♦  ‘'l) 


3u 

3y 


(130) 


3c 

Jt 
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°Wg  + 


(131) 


where 

p  —  mass  averaged  density 
u  -  mass  averaged  velocity 
i/L  -  molecular  viscosity 
c  -  sediment  volume  concentration 
W,  -  fall  velocity 
at  -  Prandtl-Schmidt  number 
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Turbulent  viscosity  is  related  to  turbulent  kinetic  energy  K  and  turbulence 
dissipation  t  by 


vt  -  Cy^e'1  0-32) 

where  is  a  constant.  The  system  of  equations  is  closed  by  the  following 

equations  for  K  and  t  : 


3K 

Jt 


(133) 


dc 

Ti 
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Sy 


♦  tC1(H  +  G) 


(134) 


where  the  indexed  values  of  C  and  a  are  experimental  coefficients,  M  is 
mechanical  energy  production,  and  G  is  buoyant  energy  production.  The 
latter  two  terms  are  equal  to 


M 


r  i 


G  ~  s  — 


(135) 


141.  Equations  130-135  are  solved  using  an  implicit  finite  difference 
scheme  using  a  time-step  of  T/120  where  T  is  the  wave  period.  Comparisons 
to  existing  data  sets  were  made  for  predicted  phase  variations  of  sediment 
concentration,  height  variation  of  period  averaged  sediment  concentration, 
maximum  horizontal  sediment  flux,  and  mean  velocity  and  shear  stress  profiles. 
Good  to  excellent  agreements  between  the  measured  and  predicted  quantities 
were  achieved  in  all  cases.  As  Hagatun  and  Eidsvik  point  out,  the  main  weak¬ 
ness  of  their  modeling  approach  is  associated  with  simulating  bed  interaction. 

142.  Glenn  and  Grant  (1987)  describe  an  eddy  diffusivity  closure  model 
that  accounts  for  suspended  sediment  stratification.  The  model  is  an  exten¬ 
sion  of  the  wave  and  current  models  developed  by  Grant  and  Madsen  (1979)  for  a 
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fixed  bed  and  by  Grant  and  Madsen  (1982)  for  movable  beds.  The  eddy  diffusiv- 
ity  closure  scheme  is  used  to  couple  the  unsteady  momentum  and  sediment  mass 
conservation  equations,  which  are  solved  for  the  combined  wave  and  current 
boundary  layer  flow,  sediment  mass  distribution,  and  suspended  sediment 
transport  rate.  The  model  also  accounts  for  ripple  formation  and  degradation, 
near-bed  sediment  transport,  and  suspended  sediment  induced  stratification. 

The  latter  process  is  represented  in  the  closure  scheme  by  applying  a 
stratification  correction  to  the  neutral,  time-invariant  eddy  viscosity.  The 
stratification  correction  term  was  found  to  be  negligible  within  the  wave 
boundary  layer.  Thus,  sediment  stratification  is  simulated  to  occur  outside 
the  wave  boundary  layer  only.  The  neutral  eddy  viscosity  is  taken  to  be 
proportional  to  the  wave-current  shear  velocity  and  varies  linearly  with 
height  above  the  bed. 

143.  Since  the  model  represents  movable  bed  effects,  and  therefore  is 
not  within  the  scope  of  this  report,  only  some  of  the  results/conclusions  from 
this  modeling  effort  are  discussed.  However,  it  appears  that  at  the  present 
time,  this  model  by  Glenn  and  Grant  is  the  most  promising  in  terms  of  repre¬ 
senting  the  pertinent  flow  and  sediment  processes  in  estuarine  boundary  layer 
flows.  Nevertheless,  the  model  is  limited  to  consideration  of  cohesionless 
sediment  transport,  and  thus  is  not  capable  of  simulating  cohesive  sediment 
transport  and  related  phenomena  (e.g.,  fluid  mud  formation/movement,  bed 
consolidation) . 

144.  Model  simulations  showed  that  the  current  velocity- induced  shear 
and  the  vertical  decay  rate  of  suspended  sediment  concentration  are  diminished 
within  the  wave  boundary  layer.  This  increased  energy  dissipation  in  the  wave 
boundary  layer  also  causes  the  current  outside  the  wave  layer  to  experience  an 
effective  (or  apparent)  bottom  roughness  that  is  larger  than  the  actual  bed 
roughness.  Above  the  wave  boundary  layer,  the  current-induced  shear  and  sedi¬ 
ment  concentration  decay  rate  are  increased,  even  more  so  by  suspended 
sediment  stratification. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


145.  These  conclusions  and  recommendations  are  interim  in  that  (a)  they 
represent  a  set  of  qualitative  opinions  based  on  the  authors'  review  of 
literature,  personal  experience,  and  judgment;  and  (b)  experimental  work  is 
ongoing  to  test  the  validity  of  the  conclusions  reached  here  so  that  a  final 
set  of  conclusions  can  be  formulated  in  a  later  report. 


Conclusions 


Estuarine  flow 

146.  Conventional  steady-flow  equations  (Equations  9-13)  for  boundary 
layer  velocity  profiles  and  bed  stresses  can  be  seriously  misleading  when 
applied  to  sediment  transport  problems  in  estuarine  flows.  They  have  been 
used  with  considerable  success  in  providing  estimates  of  water  levels  and 
general  flow  patterns  in  manual  calculations  and  two-dimensional  depth- 
integrated  numerical  models.  They  have  even  been  used  in  sediment  transport 
calculations,  but  their  use  requires  considerable  skill  and  knowledge.  Better 
descriptions  are  needed  for  detailed  (high  temporal  and  spatial  resolution) 
analysis  of  dredged  material  placed  in  open  water. 

147.  The  following  conclusions  are  drawn  from  this  review: 

a.  For  cases  where  limited  resources  prevent  the  most  rigorous 
boundary  layer  descriptions,  Equations  18-28  may  be  sufficient 
for  rapid  estimates  of  sediment  transport. 

b.  When  using  Equations  18-28,  to  take  into  account  salinity- 
induced  stratification  in  computation  of  bed  shear  (or  u*) 
and  velocity  profiles  in  estuarine  flows,  the  method  developed 
by  McCutcheon  (Equation  35)  is  recommended.  For  suspended 
sediment  stratified  flows,  it  Is  recommended  that  the  values 
for  bed  shear  obtained  using  the  method  presented  by  Coleman 
(Equation  38)  be  compared  to  those  obtained  using  McCutcheon' s 
method.  Both  methods  require  iterative  solution  of  the 
governing  equations  for  vertically  nonuniform  density  pro¬ 
files.  However,  Gust's  findings,  paragraphs  74-80,  regarding 
occurrence  of  drag  reduction  in  sediment  suspensions  (which 
causes  alterations  in  the  velocity  profile,  in  particular  in 
the  near-bed  zone)  should  be  kept  in  mind. 

c.  For  three-dimensional  modeling  and  highly  detailed  sediment 
transport  calculations,  turbulence  transport  modeling  by 
Mellor-Yamada  (Equation  37)  is  most  appropriate  and,  depending 
on  resources,  is  recommended. 
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Estuarine  flows  with  short-period  waves 

148.  The  following  methods  are  recommended  for  computation  of  bed  shear 
stresses  in  estuarine  flows  with  short-period  waves : 

a.  O'Connor  and  Yoo's  method  (Equations  65-75)  is  recommended 
when  complete  numerical  modeling  (  f  the  water  column  Is  not 
practical/possible. 

b.  Of  the  eddy  viscosity  models  reviewed,  and  considering  the 
limited  scope  of  analysis  in  this  report,  i.e.,  fixed  bed, 
Madsen  and  Wikramanayake ' s  model  (Equations  100-104  and 
113-115)  is  recommended. 

c.  The  turbulence  model  proposed  by  Davies,  Soulsby,  and  King, 
paragraphs  128-137,  is  recommended  over  other  turbulence  clo¬ 
sure  models  when  sediment  stratification  does  not  occur 
because  of  its  realistic  approach  in  defining  the  problem,  and 
because  of  its  relative  simplistic  approach  compared  with 
other  turbulence  closure  schemes . 

d.  When  stratification  does  occur,  the  turbulence  closure  model 
described  by  Glenn  and  Grant  (1987)  is  recommended. 


Re  c  ommenda  t i ons 


149.  Field  and  laboratory  data  should  be  used  to  define  the  circum¬ 
stances  under  which  these  conclusions  given  are  sound.  Numerical  experiments 
using  laboratory  and  field  data  as  guides  should  be  used  to  establish  the  com¬ 
putational  resources  required  for  the  various  levels  of  turbulence  closure 
complexxty . 

150.  After  the  appropriate  expressions  for  boundary  layer  descriptions 
over  firm  beds  have  been  selected,  they  should  be  adapted  for  soft/f lowable 
beds  using  the  general  approach  of  Mei  and  Liu  (1987). 
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Friction  Factor  Experssion  from  Theoretical  Velocity  Profiles 
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Note:  6  -  boundary  layer  thickness 

U  -  free-stream  velocity 


Table  2 


Roughness  Coefficient  Estimates 


Source  Conditions 

f 

c 

Bowden,  Fairbairn, 
and  Hughes  (1959) 

u  at  5  percent 

0.007 

Sternberg  (1968) 

u  at  100  cm 

0.006 

Tesche  (1975) 

u  at  50  cm 

0.007 

Equation  10 

n  -  0.02  ,  R  -  30  ft 

0.004 

Equation  10 

ks  -  0. 2  mm  ,  R  -  30  ft 

0 . 0004 

Equation  11 

ks  -  0.2  nun  ,  y  -  3  ft 

0.0023 

Equation  12 

u  -  3  fps  ,  y  -  3  ft 

0.0016 

Soulsby  (1983) 

Mud  bottom,  u  at  100  cm 

0.0011 

Silt  and  sand  bottom,  u  at  100  cm 

0.0008 

Plane  sand,  u  at  100  cm 

0.0013 

Rippled  sand,  u  at  100  cm 

0.0030 

Sand  and  gravel,  u  at  100  cm 

0.0012 

APPENDIX  A:  NOTATION 


Parts  II  and  III 

A  Integration  constant 

b  Channel  width 

c  Sediment  volume  concentration 

c0  Maximum  (near-bed)  sediment  volume  concentration 
ca  Depth- averaged  sediment  volume  concentration 

C  Chezy  friction  coefficient 

e  Alternating  tensor 

Ej,  Fickian  diffusion  coefficient 

f  Darcy  friction  factor 

fc  Current  friction  factor  -  f/4 

g  Acceleration  of  gravity 

G  Buoyancy  contribution  to  turbulent  energy 
h  Mean  water  depth 

i,j,k  Coordinate  axis  designator,  1,  2,  or  3  with  repeated  index 
summation 

k  Turbulent  kinetic  energy 

ka  Effective  roughness  size 

i  Mixing  length 

n  Manning  roughness  coefficient;  index  for  multiple  constituents 
in  transport;  subscript  indicating  constituent 

P„  Shear  production  of  turbulent  energy 

P  Pressure 

r  Ratio  of  mass  mixing  length  to  momentum  mixing  length 

R  Hydraulic  radius 

Rf  Flux  Richardson  number 

Rt  Richardson  number 
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5  Sediment  specific  gravity 

t  Time 

T  Period  of  oscillation 

u  Flow  velocity  in  direction  of  x  axis 

Velocity  at  specified  distance  above  the  bed 

u  Depth-averaged  flow  velocity 

u.  Shear  velocity  -  {* / f> 

u.cw  Wave-  and  current- induced  shear  velocity 
U0  Free-stream  velocity 

W,  Sediment  settling  velocity 

x  Horizontal  distance  coordinate 

y  Vertical  distance  coordinate 

y„  Bed  roughness  height 

z  Horizontal  distance  coordinate  perpendicular  to  x 

6  Boundary  layer  thickness 

Sa  Thickness  of  all  three  layers  in  Kajiura's  three-layer  model 
5l  Inner  layer  thickness 

e  Dissipation  of  turbulence 

q  Displacement  of  the  water  surface  from  mean  water  level 
8  Concentration  of  constituent 

k  Von  Karman's  turbulence  constant 

Length  scale 

v  Kinematic  viscosity  of  water 

v‘  Kinematic  viscosity  of  solid-liquid  mixture 

ut  Turbulent  eddy  viscosity 

i/tR  Reference  value  of  eddy  viscosity 

II  Wake  strength  coefficient 
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(IIA)«(y /«) 
P 
Ph 
Pti 
Pm 

p. 

Po 

a 

at 

T 


T 


sw 


* 

n 


Wake  region  velocity  augmentation  function 
Density  of  water 

Density  of  sediment-water  mixture  at  the  sediment  bed 

Density  of  sediment-water  mixture  at  the  water  surface 

Depth-averaged  density  of  sediment-water  mixture 

Sediment  density 

Reference  density 

Frequency  of  oscillation 

Turbulent  Prandtl  number 

Bed  shear  stress 

Sidewall  shear  stress 

Sum  of  sinks  and  source  for  9 

Coriolis  effect  parameter 


Superscripts 

'  A  quantity  fluctuating  in  time 


Subscripts  to 

r  and 

S 

xy 

In  the 

X 

direction  on  a  y 

constant  surface 

xz 

In  the 

x 

direction  on  a  z 

constant  surface 

Instantaneous  value 
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Part  IV 


a  Amplitude  of  water-surface  motion 

Ao  Wave  orbit  horizontal  amplitude 

A*,*,  Wave  orbit  horizontal  amplitude  at  the  bed 
c  Sediment  volume  concentration 

C  Chezy  friction  coefficient 

d  Height  of  inner  layer  of  two-layer  viscosity  model 
fc  Current  friction  factor  -  f/4 

f„  Wave  friction  factor 

fwc  Friction  factor  for  combined  waves  and  currents 
g  Acceleration  of  gravity 

G  Buoyant  turbulent  energy  production 

h  Depth  of  flow 

k  Wave  number  -  2w/L 

kb  Characteristic  dimension  of  bottom  roughness 

k„  Boundary  roughness  size 

ka  Equivalent  bed  roughness  size 

K  Turbulent  kinetic  energy 

L  Wave  length 

£  Mixing  length 

£'  Measure  of  wave  boundary  layer  thickness 

M  Mechanical  energy  production 

p  Pressure 

vp  Mean  pressure  gradient 

t  Time 

T  Period  of  oscillation 
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Tyi  Reynolds  shear  stress  in  the  z-direction  associated  with  the 
current 

u  Velocity  in  direction  of  x  axis 
u  Mass  averaged  velocity 

ux  Velocity  in  x-direction  at  specified  distance  above  the  bed 

u*(t)  Magnitude  of  modified  shear  velocity 

b 

ubm  Maximum  wave  velocity  at  the  top  of  the  wave  boundary  layer 
ud  Wave-induced  orbital  velocity  within  wave  boundary  layer 
u„  Wave  orbital  horizontal  velocity 
u^  Maximum  wave  orbital  velocity  near  the  bed 
u,.  Resultant  velocity  of  simultaneous  wave  and  flow 

u.  Shear  velocity  -  yr/P 

11*,  Near  bottom  velocity 

uuc  Shear  velocity  associated  with  the  current  boundary  shear  stress 

|u.c|  Magnitude  of  the  time-averaged  shear  velocity 

|u,wc|  Magnitude  of  maximum  shear  velocity  due  to  combined  wave  and 

current  motion 

U«  Amplitude  of  the  horizontal  wave-induced  velocity  component 
immediately  above  the  wave  boundary  layer 

U  Depth-averaged  current  speed 

|Ua  |  Magnitude  of  steady  current  velocity  at  height  y  -  a  above  the 
bottom 


U0  Free-stream  velocity 
0  Laminar  current  velocity  at  y  -  5/(tr/2v)0-5 

v  Velocity  in  direction  of  y  axis 

w  Velocity  in  direction  of  z  axis 

Wg  Sediment  settling  velocity 
x  Horizontal  distance  coordinate 
y  Vertical  distance  coordinate 
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y'  Hypothetical  elevation  above  the  bed  (proposed  by  Bijker)  where 
a  line  emanating  from  the  seabed  (y  -  h)  is  tangential  to  the 
current-only  velocity  profile 

yx  Height  above  the  theoretical  bed  at  which  the  velocity  is  equal 
to  zero 

y„  Bed  roughness  height 
z  Horizontal  distance  coordinate 
a  Current  velocity  reduction  factor 
a'  Ratio  of  eddy  diffusivities  of  turbulent  energy  and  momentum 
6  Boundary  layer  thickness 
6W  Wave  boundary  layer  thickness 

«  Turbulence  dissipation  rate 

r)  Displacement  of  the  water  surface  from  mean  water  level 

9  Angle  between  orbital  wave  motion  and  the  direction  normal  to 
the  current 

k  Von  Karman's  turbulence  constant 

H  Relative  velocity  of  maximum  wave  orbital  velocity  to  the 
current  at  y' 

v  Kinematic  viscosity  of  water 

Molecular  viscosity 

i/t  Turbulent  eddy  viscosity 

p  Density  of  water 

p  Mass  averaged  density 

a  Wave  angular  frequency 

ac  Turbulent  Prandtl  number 

r  Bed  shear  stress 

rb  Bed  shear 

rbm  Maximum  wave-  and  current-induced  bed  shear 
rc  Bed  shear  stress  due  to  current  only 

r.  Maximum  combined  wave-current  bed  shear  stress  vector 
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r„  Maximum  wave  bottom  shear  stress 
rwc  Resultant  combined  shear  due  to  waves  and  currents 
rb  Current-induced  time  mean  bed  shear  stress 
rwc  Time  mean  wave  and  current  induced  bed  shear  stress 
<(>  Angle  between  current  and  wave  motion  directions 
4>c  Angle  of  the  average  shear  stress 
Superscripts  to  r  and  6 

’  Surface  (grain  or  net)  shear 
"  Form  shear 

No  superscript  indicates  total  shear 


Subscripts  to 

r  and  6 

c 

Caused  by 

current 

w 

Caused  by 

waves 

wc 

Caused  by 

waves  plus  currents 

xy 

In  the  x 

direction  on  a  y  constant  surface 

xyo 

In  the  x 

direction  at  the  bed  (y  -  0) 

xz 

In  the  x 

direction  on  a  z  constant  surface 
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